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Abstract 



Using an error models motivated by the Knill, Laflamme, Milburn proposal for efficient linear 
optics quantum computing [Nature 409,46-52, 2001], error rate thresholds for erasure errors 
caused by imperfect photon detectors using a 7 qubit code are derived and verified through simu- 
lation. A novel method - based on a Markov chain description of the erasure correction procedure 
- is developed and used to calculate the recursion relation describing the error rate at different 
encoding levels from which the threshold is derived, matching threshold predictions by Knill, 
Laflamme and Milburn [quant-ph/0006120, 2000]. In particular, the erasure threshold for gate 
failure rate in the same order as the measurement failure rate is found to be above 1 .78%. 
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Chapter 1 
Introduction 



Quantum computation was bom when Benioff Q proposed using the laws of quantum mechanics 
instead of classical mechanics to perform computation. More importantly for physicists, Feyn- 
man |l6| pointed out, was the possibility of simulating quantum system by using such quantum 
computers - because the state space of quantum systems has dimensions that grow exponen- 
tially with the number of subsystems involved, they are notoriously inefficient to simulate with 
classical computers. 

Much progress has been done on how to perform certain algorithms much more efficiently in 
quantum computers, as well as how to simulate quantum systems with quantum computers. The 
main challenge is to construct physical systems on which such computers may be built. Many 
proposals have been put forth, each with its strengths and weaknesses, but the recent proposal 
by Knill, Laflamme and Milburn ifTSl . which uses only linear optics elements and single photon 
sources and detectors, is of particular interest for various reasons. First of all, it was believed that 
it was impossible to build a universal quantum computer from linear optics elements because 
photons do not interact with each other. This proposal relies heavily on state preparation in order 
to avoid such a hurdle. Moreover, one of the most successful applications of quantum comput- 
ing, quantum key distribution, would benefit directly from the construction of an optical quantum 
computer because most protocols are implemented with optical communication. Quantum com- 
putation with non-linear optics, on the other hand, is notorious for the photon loss, making it 
extremely inefficient. 

The main problem in a physical implementation of a quantum computer is its sensitivity to 
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error and to the environment. In a classical computer, information is represented by two different 
states, and i, and such computers can easily be built so that error in distinguishing between 
these two states is insignificant. In a quantum computer, the state may be in a superposition 
of these two basis states, say (x|o) + |3 and when the state is measured, we obtain |o) with 
probability |ap, and |i) with probability |(3p. The information in a quantum computer ends 
up being more akin to an analog data than to digital data, even though there is only a discrete 
number of basis states. Imprecision in devices used for computation, the logic gates, has a much 
greater impact in quantum information than in classical information. Moreover, it is much harder 
to isolate a quantum system from the environment, so that interaction between the two also end 
up introducing error into the computation, in an effect known as decoherence, which is only 
observed in quantum systems. 

The linear optics quantum computing proposal also suffers from these problems, even though 
photons are much less sensitive to decoherence than other physical systems proposed for quan- 
tum computation. The reason is that the gates proposed for efficient linear optics quantum com- 
putation are probabilistic, so that even with infinite precision in the linear optics elements, they 
may fail. The authors of the proposal have demonstrated that by encoding the data carefully, 
one can easily overcome the probabilistic nature of the gates ifTTI . However, these gates depend 
on single photon measurements, which are known to have low efficiency, although the gates 
described in the proposal can flag when detectors have failed and replace the lost qubit automat- 
ically. The natural question then is what is the maximum rate at which photons may be lost in 
these gates while still allowing for scalable quantum computation? The main objective of this 
thesis is to answer that question. 

This thesis is organized as follows. Chapter 2 gives a basic overview of quantum erasure cor- 
rection codes and fault-tolerant computation, focusing on CSS codes and the stabilizer formal- 
ism. Chapter 3 describes the basics of the efficient linear optics quantum computation proposal 
by Knill, Laflamme and Milbum, focusing on the description of the error model for teleportation 
failures due to the inherent probabilistic nature of the linear optics implementation, as well as 
photon loss due to detector inefficiencies. Chapter 4 gives a detailed description of the Steane 
code, and contrasts it to the Grassl code in the context of the error model for linear optics quan- 
tum computing. Given these background chapters, it becomes clear how to perform the erasure 
correction procedure, and in Chapter 5 a new and compact description of the procedure is given 
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in terms of Markov chains. With this description the threshold is obtained for both error models, 
with the detailed calculations being presented in the appendices. In Chapter 6 a Monte Carlo 
simulation of the erasure correction procedure is described, and the simulation results are pre- 
sented and contrasted with the theoretical prediction of Chapter 5. The thesis concludes with a 
discussion of the results, and suggestions for future work. 



Some basic notation is assumed in this thesis, and it is briefly reviewed here. Other notation is 
introduced in the body of the thesis as needed, along with explanations. 

In order to emphasize the difference between qubits and photon number states, qubits in the 
computational basis will be represented with gothic font, e.g. |o) and while photon numbers 
will be represented in the usual fonts, e.g. 11) and |2). Operators acting on the physical qubits 
(qubits not encoded by error correction code) are denoted by bold capital letters, e.g. X, and 
multi-letter operators use typewriter font, e.g. CNOT. Operators acting on the encoded qubits 
are denoted in the same way, with an added overbar, e.g. X, and encoded qubits also have an 
overbar, e.g. |¥) or |oT). If it is clear from the context that the operator is encoded, the overbar 



1.1 Notation 



is omitted. 



The single qubit Pauli group =V consists of the operators 



Y 



X 



Z 



I 




(1.4) 



(1.2) 



(1.3) 



(1.1) 



along with multiplication by i = v^— T, following the convention in ifTOl . Under this definition. 
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products of Pauli operators are given by 



xz 



(1.5) 
(1.6) 
(1.7) 
(1.8) 



The Pauli group Vn over n qubits is given by the n fold tensor product of single qubit Pauli 
operators along with multiplication by {±i}. 



Other common unitary operations are 
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(1.9) 
(1.10) 

(1.11) 



(1.12) 



where H is the Hadamard transform, P is the phase gate, CNOTi 2 is the controUed-X where 
qubit 1, the control, determines the application of an X on qubit 2, the target, and CSIGNj 2 is 
the controUed-Z with control and target similarly defined. Throughout the thesis, H and P are 
taken to be a single qubit operations. 

These unitary operations can also be represented by rotations about Pauli operators. In gen- 
eral, given some operator U G Vn, we define 



(U) 



(1.13) 
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where 9 is given in degrees. In particular, if 9 = ±180°, then the rotations are equivalent to 
operators in Vn, and if 9 = ±90° the rotation is equivalent to a product of H, P and CNOT 
applied to various different qubits. 



Chapter 2 

Erasure Correction and Fault- Tolerance 



A brief overview of erasure errors is given, along with a brief discussion on the basics of quantum 
error correction codes and fault-tolerant computation. In this discussion, quantum error correc- 
tion codes are referred to as error correction codes. Error correction codes over classical data 
will be referred to as classical error correction codes. 

2.1 The Erasure Channel 

A common error model for quantum data is given by a channel for which there is a finite proba- 
bility Pr(0) of an error superoperator O being applied to the qubit transmitted. In that case, we 
can define the channel by the superoperator 

£M = (1 -Pr(0))p + Pr(0)0(p), (2.1) 

where p is the density matrix of the qubit input into the channel, and we take this channel to be 
memoryless, that is, different uses of the same channel are statistically independent. 

In general the corruption of data is not a priori obvious to the observer, and as was described 
in the introduction, one must encode the data in special ways in order to detect such corruption. 
Under some physical models, however, it is immediately known when some error superoperators 
have been applied. The canonical example of this is spontaneous emission in qubits represented 
by atoms lfT2l . where one may detect the resulting photons and determine that the state of the 
atom has been corrupted. In general, it is possible to consider the qubit encoding to be a Hilbert 
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space strictly smaller than the Hilbert space describing the entire physical system. Say for exam- 
ple, the computational Hilbert space is given by 

= span{|o) , (2.2) 

but the state of the entire physical system is in the Hilbert space 

n = span{|o) , |i) , |2) , • ■ ■ , In)}. (2.3) 

Errors that map qubits into the space orthogonal to Hq, in a process known as leakage, can 
always be detected without disturbing the computational subspace, and therefore can be consid- 
ered erasures, as long as the qubits determined to be in Hj^ are replaced by fresh computational 
qubits. This is, in essence, the case for the linear optics proposal described in Chapter^ and it 
is the main motivation for this thesis. 

Abstracting from the implementation details, we can think of an erasure channel as an er- 
ror channel with side information through which the application of some error superoperator is 
flagged, and therefore it is known which qubits have been corrupted. In that case, because there 
is certainty that the error superoperator has been applied, we call it an erasure superoperator, 
and we represent it by a hollowed letter corresponding to the error superoperator applied, or O. 
The important distinction being that it is known with certainty that the operator has been applied, 
so the probability of the error O occurring is absent from this description. 

The superoperator equivalent to a full qubit erasure, or a complete randomization of the state, 

is 

E(p) =I = l(p + XpX + YpY + ZpZ]. (2.4) 

The partial erasure superoperators, which correspond to partial randomizations of the qubit (a 
unique quantum feature) are given by 

Z(p)=l(p + ZpZ) (2.5a) 

X(p]=l(p + XpX), (2.5b) 

corresponding to a possible X error or a possible Z error respectively, and called X erasure 
and Z, or phase, erasure. As will be demonstrated in Chapter |3l photon loss due to detector 
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inefficiencies within a CSIGNs in linear optics quantum computing can be modeled as a gate 
application followed by transmission over an erasure channel. 

For the rest of this thesis, we will refer to erasure superoperators that act independently on 
multiple qubits as erasure patterns. For example, a block of seven qubits where the first is 
affected by a Z erasure, while the last is affected by a full erasure, is said to have been affected 
by the erasure pattern 

Z(g)I(g)I®I®I®I®E. (2.6) 

This, in effect, describes a convex sum' of Pauli operators acting on seven qubits. These Pauli 
operators are known as error operators, and they represent the actual errors that are detected by 
the error correcting procedure. In our discussion, we will call them erasure operators when it 
is known which qubits could have been corrupted. In particular, given (12. 6I ). measurements may 
indicate that any erasure operator among 
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(2.7) 



was applied, and given the types of erasures considered here, the probability distribution for any 
of these erasure operators being measured is always uniform. By definition, the weight of an 
n qubit Pauli operator is the number of qubits over which it acts non-trivially. Similarly, the 
number of qubits over which some erasure pattern W may act non-trivially is called the weight 
of the pattern, denoted wt(W]. 

'in this context, we take a convex sum to mean the application of different Pauli operators Oi with probability 
Pr(Oi) > to a state described by some density matrix p so that Y.i Pr(Oi) = 1 , yielding some superoperator 

0(p) =^Pr(Ot)OtpO|. 
It is clear that ^2A\ and ^2.5\ fit into this category. 
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2.2 Conditions for Quantum Erasure Codes 

There are special conditions for an encoding of quantum data to be considered a quantum error 
correcting code, and they hinge on what error operators a code claims to be able to correct. 
The conditions developed by Knill and Laflamme lfT6l are, given a quantum code consisting of 
encoded states {|ct)}|^^^ and correcting a set of error operators {Ej} 

(C||E|Ej|ct) = (cJEjEjlc^) (2.8a) 
(c;|E|Ej|c^) = for (ctlc^) = 0, (2.8b) 

which are known to be necessary and sufficient. If the E| are Pauli operators with maximum 
weight t, then this code is a distance 2t + 1 code. 

The knowledge of exactly which qubits have been corrupted is very powerful, and in general 
it allows for twice as many single qubits to be corrupted while still allowing the data to be 
recovered. Therefore, given some general error correcting code with parameters^ [[n, k, d = 
2t + 1 ]], one will only be able to correct t errors at unknown locations, while being able to correct 
at least It erasures^. This follows from the fact that there is no need to distinguish between error 
operators at different locations because the positions of the corrupted qubits are known when 
dealing with erasures. With that in mind, Grassl et al. lfT2l derived modified conditions for 
quantum erasure codes. In the general Knill-Laflamme conditions, each of the error operators 
are taken to have weight up to t, so that the product E|E^ has weight up to 2t. Condition (I2.8al) 
therefore says that valid states disturbed by an error operator E|Ej are scaled and rotated in the 
same way. Condition (I2.8bl) on the other hand guarantees that E|Ej will never map orthogonal 
encoded states into non-orthogonal states, so that the different basis states, when corrupted, are 
still perfectly distinguishable. Since it is known over which qubits E|Ej acts non-trivially, this 
erasure operator can also be corrected. Thus one can talk about a set of correctable erasure 
operators made up of all products {E|Ej} where E,,, Ej are correctable errors. If the set of erasure 
operators is relabeled {At) then they satisfy the Modified Knill-Lafiamme Conditions lfT2l 

(ctlAtlct) = (cJAilc^) (2.9a) 
(ctlAJc^) = Ofor (ctlc^) =0, (2.9b) 

^That is, it encodes k qubits using n qubits and being able to correct Pauli errors of maximum weight t. 
■'The same can be said about classical error correcting codes. 
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making it clear that t error correcting codes are in fact 2t erasure correcting codes. Note that 
the same general method is employed to identify which underlying Pauli operator acted on the 
data regardless of whether it is an erasure or a general error, as will be demonstrated later in 
this chapter. For an erasure pattern, however, the qubits over which an operator may have acted 
non-trivially are known, so all the are only required to come from the same erasure pattern. In 
essence, independently for each correctable erasure pattern, all the Pauli operators in the convex 
sum describing the erasure pattern must satisfy ( I2.9al) and (I2.9bl ). 



One of the first large classes of quantum codes to be discovered were codes based on pairs 
of orthogonal classical codes. These codes are named CSS codes in honor of the discoverers: 
Calderbank, Shor 123 and Steane Q. 

To understand how CSS codes are constructed, first we need to quickly review some basic 
facts about classical error correcting codes ll26ll . 

A classical code C is a set of length n vectors over GF(2) (the integers modulo 2). Each of 
these vectors is called a codeword, and if this set forms a subspace of GF(2)"^ then it is called 
a linear classical code. From here on classical linear codes will be referred to simply as linear 
codes. The Hamming weight 'wt(c) of a vector v is the number of non-zero elements. The 
Hamming distance of two vectors x,y is defined by wt(x — y). Recall that, because we are in 
GF(2), 1+1 =1—1 = 0, so the operation can be thought of as an elementwise XOR. 

A linear code C with 2^ codewords of length n is said to be a [n, k, d] code if the minimum 
weight of a non-zero codeword is d - this is also the minimum distance between two distinct 
codewords. Because this is a linear subspace of dimension k, we can describe the code by a 
generator matrix G with dimensions k x n, where the rows are linearly independent codewords, 
and then every uncoded length k row vector v can be encoded into a codeword c by performing 
the matrix multiplication 



The same code C may be described implicitly by the codespace orthogonal to it, C . That is, we 
can construct a n — k by n matrix H such that 



2.3 Calderbank-Shor-Steane Codes 



vG 



= c. 



(2.10) 



He = 



(2.11) 
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if c is any valid codeword in C. 

The minimum distance of a linear code is important because any vector over GF[2)"^ can be 
associated with at most one valid codeword if it is within a Hamming distance of 

d- 1 



2 



(2.12) 



and therefore we call a code with distance d a t error correcting code. If all we are interested in 
is detecting whether an error has occurred or not, we can tolerate at most d — 1 errors, since d 
errors may lead one valid codeword into a different one without being detected. 

Since it is known that He — for any valid codeword, if we consider some error vector e, 
then we have 

H(c + e] =Hc + He = He. (2.13) 

As long as e ^ C, this is a non-zero value, and it will tell us that an error has occurred - we have 
detected an error. If we assume 'wt(e) < t, then the value He is uniquely mapped to e, and one 
may simple apply e to the corrupted data to obtain c + e + e = c. 

The n — k rows of H are linearly independent vectors that are orthogonal to all valid code- 
words in C, as previously stated. So these vectors may be thought of as a basis for the orthogonal 
codespace C-*- (which may share more than the zero codeword with C). 

Consider the case where C"*" C C are linear codes with n bit long codewords'*. We say that 
two codewords in C are equivalent if they differ by an element of C^. If C is a [n, k, d] code, 
then encodes n — k bits, and by a simple counting argument, it is clear that there must be 
2k — n equivalence classes of C-*- in C. We define the quantum superposition |ic), for some 
c e C, to be 

\^c) = ^= }2 Ic + u). (2.14) 

This is, in effect, a superposition over the equivalence class containing c. It is clear that if Ci and 
C2 are non-equivalent classical codewords 

(icjic,) -0, (2.15) 

and since there are 2k— n of these, we consider the space spanned by collection of all possible |ic) 
to be a [[n, 2k — n]] quantum code. Since the codewords of C are guaranteed to have minimum 
^This is not the most general CSS code construction, but for simplicity we restrict ourselves to this case. 
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distance d, |ic) will always be orthogonal to a state resulting from the application of less than d 
bitflips, or X operators. 

Applying the qubitwise Hadamard transform to (12.141) one obtains^ 

|tc> = -|=}l(-1)""|u), (2.17) 

which is a superposition over codewords of C with relative phases dependent on c. Again, 
because C has minimum distance d, if we apply less than d bitflips to |ic) we obtain a state 
orthogonal to it, which is not the case if we apply d bitflips. Because we are in the Hadamard 
basis, that means that the |ic) basis can withstand at most d — 1 phase flips, or Z operators, and 
the error will still be detected. By analogy, we define this to be a [[n, 2k — n, d]] quantum code, 
since there are weight d Pauli operators that cannot be detected as errors. 



2.4 Stabilizer Codes 

Developed by Gottesman|l8l|, and by Sloane, Shor, Calderbank and Rains ll25ll under a different 
formalism, stabilizer codes are a class of quantum codes much broader than the one described by 
CSS codes. 

The basic idea is to define an Abelian subgroup S of the n qubit Pauli group Vn, and to take 
the common eigenspace with eigenvalue +1 of S as the code space. S is called the stabilizer of 
the code, since any operator M e S and encoded state gives 

M|{b> = |M^>, (2.18) 

that is, all elements of S act trivially on the codespace. It can be shown that if S has n — k 
generators, then the codespace has dimension 2^, and therefore it encodes k qubits. 

Now consider an error operator E G Vn- If for some M G S we have^ {M, E} = 0, then 

ME = -EM = -E jijJ) (2.19) 



By using the identity 

ifu^C 



CGC 

^Recall that {A, B), the anti-commutator between A and B, is defined as {A, B} = AB + BA. 
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so E lij)) is a — 1 eigenstate of M, and this can be detected by measuring M. Going back to the 
modified Knill-Laflamme conditions for erasure correction, one finds that if At G S, then (I2.9al ) 
holds, since 

(ct|Adci) = (cdcO = l, (2.20) 
for any encoded states |ci,). If (At, M} = for some M G S then (I2.9bl) holds, since 

{ci\ Ai |c^) - (cil A|M |c^) = - {ci\ MA| |c^) = 0. (2.21) 

so we can restrict ourselves to looking at the commutation relations of erasure operators to de- 
termine their correctability. 

2.4.1 CSS Codes as Stabilizer Codes 

From the basis state representation of CSS codes, we can infer what the stabilizer of the code 
should be. Recall that the equivalence relation that defines the |tc) in (12.141 1 is that if two code- 
words in C differ by an element of C"*", then they are equivalent. Clearly, if we add some w G C"*" 
to each of the classical codewords in |ic) the state remains unchanged - all codewords remain 
in the same equivalence class. We can take this operation to be a bit flip operator Wx, obtained 
by replacing each in w by an I, and each 1 is replaced by an X, all elements concatenated by 
tensor products. This is a set of bit flip operators in the stabilizer, and they can be generated by 
the n — k linearly independent bitflip operators obtained from the n — k generators of C-*- . 

Looking at the Hadamard basis description of the CSS codes as given in (12.171) . we notice a 
similar fact, with a little more algebra involved. Again, if we apply a bit flip operator Wx based 
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on a codeword w G C"*-, we obtain 



Wx|ic> = ^}2(-1)^-"|u) (2.22a) 



CI 

1 



itec 



CI 

1 



CI 

1 



CI 

1 



^(-1)^-"!u + w) 

itec 

u-'ec 

^(_l)C.u'-c.W|^/^ 

u'ec 

^i-^r^' iu') 

it'ec 

Wx|tc> = |tc> (2.22b) 

Thus Wx stabilizes the |ic) as well. In the computational basis, this operator is the same, except 
that the Xs are replaced by Zs, so we write it as the operator Wz. Again, there are n — k linearly 
independent operators such as these, since there are n — k generator codewords for C"*". 

Taking all the possible Wz and Wx, we have 2n — 2k linearly independent generators, and 
we have n qubits. Thus, with these generators we encode n — [In — 2k) = 2k — n qubits, 
which is exactly the number of qubits that the CSS code encodes, so the Wz and Wx generate 
all stabilizer operators of the CSS code. 

For any self-orthogonal linear code C with parity check matrix H, we obtain a CSS code with 
stabilizer generators obtained from H as described in such a way that we have n — k generators 
that are tensor products of only Xs and Is (what we call the X stabilizers), and n — k generators 
that are tensor products of Zs and Is (what we call the Z stabilizers). 

2.4.2 The Normalizer and the Heisenberg Representation 

There are Pauli operators that commute with all elements of the stabilizer S but that do not 
necessarily leave the codespace invariant. That is, there are operators O G Pn such that 

VM G S, [M,0] = (2.23a) 
0|^) = (x|x), (2.23b) 
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where |\[>) , |x) are encoded basis states (possibly equal) and a = ±1 . The set of such operators 
is called the normalizer of S, denoted N(S). Operators in N(S)/S are errors that cannot be 
detected because they map valid codewords into different valid codewords, and because of that 
they can also be seen as encoded operations on the encoded data - in fact, they are the encoded 
Pauli operations acting on encoded qubits. 

Observing the evolution of S and N(S] under the action of different unitary operators can 
be used to determine the behaviors of certain types of circuits, and it is especially helpful in 
constructing encoded operations lfT0ll33l . This is what is called the Heisenberg representation 
of quantum computers |l9|, since it is based on the general idea of tracking the evolution of 
the operators in N (S) - S simply being a particular subset of N (S) - much like one tracks the 
evolution of operators in the Heisenberg picture of quantum mechanics. The general idea is to 
observe how some operators M| G N(S) evolves under the action of some unitary operator U, 
by noting in particular that 

UMji|;) = UM|U"^U (2.24a) 
UMtMjU"^ = UMiU"fUMjU^ (2.24b) 

We may interpret (I2.24all as a statement of what Mt is mapped to under the action of U, so that, 
for example, we may know how its eigenstates get mapped under the action of U. On the other 
hand, (I2.24bl) allows us to restrict our attention to a generating set of unitary operations - the 
action of the generated group follows by linearity. 

However, in general, a unitary operator will not map a Pauli operator into another Pauli 
operator. We can consider, however, a particular class of unitary operators that is very useful. 

Definition 1. The Clifford Group, denoted 62, is the set of unitary operators that maps the Pauli 
group into itself under conjugation im\l . That is 

62 = {UIUOU^ G Vryfor all O G Vn}. (2.25) 

Since S C Vn and N(S) C V-n, we can consider circuits made up of gates in 62, and the 
Heisenberg representation allows us to monitor the evolution of the encoded states by observing 
the evolution of operators that generate the encoded Pauli set "P^. 

C2 is finitely generated by the Pauli group plus the Hadamard, H, the phase gate P = 
diag(l,i), and the CNOT. If we restrict our attention to the action of these gates, we easily 



16 



Erasure Thresholds for Linear Optics 



find how any normalizer is transformed. For the purposes of this thesis, it suffices to look at the 
action of the generating set of unitary gates of 62- 
The Hadamard gate maps the Pauli operators as 

HXH = Z (2.26a) 
HZH = X, (2.26b) 

so that, for example, we may consider an X followed by a H to be the same as a H followed by a 
Z. This, in fact, is a very useful tool in observing the propagation of error operators in quantum 
circuits. The mappings for the other gates in the generating set of 62 are 

PXPt = lY (2.27a) 

PZP"^ = Z (2.27b) 

CN0Ti^2(X® I)CN0Ti,2 = X(g)X (2.27c) 

CN0Ti_2(I® X)CN0Ti_2 = I®X (2.27d) 

CN0Ti,2(Z ® I)CN0Ti,2 = Z®I (2.27e) 

CN0Ti,2(I®Z)CN0Ti,2 = Z(g)Z, (2.27f) 

along with the identity 

(I(g)H)CN0Ti_2(I®H) = CSIGNi_2, (2.28) 

will be used extensively in the construction of fault-tolerant encoded gates in Chapter |4l Note 
that if we consider instead rotations about Pauli operators, then all 90° rotations and its integral 
powers will also be in Gj- These, however, are just a different representation of the gates that can 
be generated by the set of gates described above. 



2.5 Erasure Correction Procedure 

The process of erasure correction hinges on the measurement of the stabilizer operators for the 
code, like in error correction codes, but now there is the extra knowledge of where the erasure 
occurred. Here we follow an approach proposed by Zalka ll32ll and based on the work of Shor 
lE^ . The general idea is that, given a single erasure in an erasure pattern which may include 
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several, one attempts to measure a stabilizer that acts trivially on all other erasures, but that acts 
non-trivially on the erasure that is being targeted for correction - what is called the correction 
target - and on qubits unaffected by erasures. The outcome of the measurement, which is de- 
scribed in more detail in the next section, indicates which action must be taken to correct the 
erasure. 

In the case of a Z erasure, it is necessary to determine whether a Z error has indeed been 
applied, or if the identity has been applied. Thus, it is sufficient to measure a stabilizer operator 
with an X on the same position as the correction target. In the case of a full-erasure, the erasure 
needs to be corrected in two separate steps. The reason for that is clear by noting that 

X(Z(p)) = l[Z(p)+XZ(p)X] (2.29a) 

= l(p + XpX + ZpZ + YpY] (2.29b) 
= E(p) (2.29c) 
= Z(X(p)), (2.29d) 

so we may consider a full erasure to be an X erasure followed by a Z erasure. First we measure 
a stabilizer with a Z on the correction target position, and correct for an X erasure, and then we 
measure a stabilizer with a X on the correction target position, correcting for a Z erasure in that 
position. 

When considering a single erasure, there is no difference between choosing to correct the 
Z erasure or the X erasure first. However, when considering an erasure pattern that consists of 
multiple erasure of multiple type, the best strategy is to make a choice that is least likely to lead 
to an uncorrectable error. This depends both on the error model and on the code being used, and 
will be discussed in more detail in Chapter|4l after both topics have been introduced. 



2.6 Fault-tolerant Computation 

We would like to be able to perform useful computation on a quantum computer regardless of 
how long the computation is or how many qubits are involved, simply because we would like to 
solve many different types of problems, of different complexities, with different input sizes. If 
one expects the error rate of the quantum computer to be naturally low enough so that errors are 
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unlikely to occur during computation, one finds that the acceptable error rates are dependent on 
the size of the computation. Thus, we'd like a means to perform any useful quantum computation 
even in the presence of a fixed probability of error for each gate. This is what is generally meant 
hy fault-tolerant quantum computation. 

Encoding the data to resist error operators is not enough to reach this objective. If, in order to 
do computation, one was required to decode the data, perform computation, and then re-encode 
the data, there would be no protection from the noise and decoherence during the computational 
step. Universal computation with the protection of error correction codes is not trivial, however, 
since it requires that the encoded operations necessary for universal computation be identified. It 
is not even sufficient to perform these encoded operations correctly, because it is possible that the 
computation still allows errors to propagate in a catastrophic way - it could be that even during 
a step where no errors have occurred, careless computation could take a correctable error into an 
uncorrectable error. The canonical example of this is the CNOT gate. It has been stated before 
that CNOT is in the Clifford group, and the exact mapping it performs - as seen in the previous 
section - does not preserve the weight of length two Pauli operators. Clearly, if we perform the 
CNOT between two qubits of the same code block, we run the risk of increasing the weight of 
the error, possibly leading to an uncorrectable error, even when none of the CNOTs fail, simply 
because the data contained errors that were propagated carelessly. A general rule that can be 
extracted from this is that we do not want errors to propagate within a code block, so we do not 
allow for qubits in the same code block to interact with each other. This, essentially, translates 
to the requirement that encoded gate operations be transversal - that they operate qubitwise on a 
code block. 

2.6.1 Fault-tolerant Stabilizer Measurement 

In stabilizer codes, in order to determine which error has affected the data, one needs to measure 
some subset of stabilizer operators, and in general the stabilizer generators suffice. According 
to (12.191 ). we know that a detectable error has eigenvalue —1 . On the other hand, the absence of 
errors or undetectable errors have an eigenvalue +1 , so we can use the phase kick-back quantum 
circuits to measure the eigenvalue of the data, as depicted in Figure 12.11 The control qubit will 
be in the state |o) + |i) after the Hadamard gate. Because M is a stabilizer of the codespace, 
whether it is applied or not does not affect valid data at all, and if there is no error (or if the 
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Figure 2.1: Quantum circuit for measuring the eigenvalue of a stabilizer operator. 



Ml 



M, 



M2 



Figure 2.2: Quantum circuit for measuring the eigenvalue of a weight three stabilizer operator 
M = Ml ® M2 ® M3 fault-tolerantly. This generalizes in the obvious way for stabilizers of 
higher weight. 



error commutes with M), both control and data are unaffected, the second Hadamard brings the 
control back to |o) and that is the state that is measured. If there is an error that anticommutes 
with M, then the control will be in the state |o) — |i), and the second Hadamard will bring it to |i) 
which is then detected. Thus, clearly detecting a |o) indicates commutation, while a |i) indicates 
anticommutation . 

Because the data is made up of multiple qubits, and the stabilizer acts non-trivially over more 
than one of them, the circuit in Figure 12.11 is not fault-tolerant. This is because if one of the 
CNOTs or CSIGNs (for Xs and Zs in M respectively, not shown in the figure) fails, it can cause 
an error on the control bit, which is shared by all the other controlled operations on other qubits. 
For example, if a failure occurs on the first controlled operation, all other qubits over which M 
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acts non-trivially will be affected. 

The solution for this problem is quite simple. If we replace the single control qubit by a group 
of qubits, one for each controlled operation needed to measure the stabilizer M, and interact with 
the data transversally as depicted in Figure IT21 there is no catastrophic error propagation. In this 
case, the control lines need to be replaced by a 'cat' state, of the form 

|o)®™+|i)®™, (2.30) 

where m is the number of qubits over which M acts non-trivially. If the error anticommutes with 
M, this cat state will develop a relative phase of — 1 . 

\of^-\if^. (2.31) 

Recall, however, equation (12.171) . and take C to be the ra bit repetition code - that is, a code that 
maps 

^ 00^ (2.32a) 

m times 

1 ^ IL;;_1- (2.32b) 

m times 

If we apply the qubitwise Hadamard to (I2.31I) . we obtain a superposition of odd weight binary 
strings, and by measuring and computing the parity classically, we can fault-tolerantly detect that 
an error has occurred. 

The use of cat states is necessary because it ensures that no information about the state of 
valid data is transfered to the ancilla. This is important because by measuring the ancilla we do 
not want to collapse the superposition of valid encoded states with no errors. 

2.7 Threshold Theorem 

Once the data is protected by an erasure correction code, the probability of erasure on the en- 
coded data may still be unacceptably high. One way to get around this problem is to perform 
concatenated coding, that is, encode the encoded data once again. Say, for example, we have 
a code C*^' with parameters [[n, 1 , d]]. By concatenating the code with itself once, we obtain 
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a [[n^, 1 , d^]] code which we call C*^', on which erasure correction can be seen as erasure cor- 
rection on each level of concatenation separately. It is straightforward to see why the minimum 
distance scales in such a way. If C has as its encoded operators X and Z, then by replacing the 
uncoded Pauli operators by these encoded ones in the stabilizer of C, we obtain a new code C*^', 
which also includes the n-fold Cartesian product of the stabilizer of C. This procedure can be 
repeated L times to obtain the concatenated code C''^' with parameters [[n^, 1 , d^]]. 

Consider C'^' once again. If the qubit erasure rate is e, then we can write the block failure 
rate as 



N 



^Cie\ (2.33) 



i=d 

where Ci are integer coefficients dependent only on the erasure correction procedure and i. We 
are only interested in the case d > 2, since we want to be able to correct at least one erasure. This 
recursion relation can be calculated by analyzing how erasures are introduced during the erasure 
correction procedure, and how this may lead to an uncorrectable erasure pattern in a code block, 
what we call a failure. In theory, one often assumes that erasure correction is attempted until 
the data is erasure free, leading N to be infinite. In practice, however, only a certain maximal 
number of erasure correction steps are attempted, placing a bound on N . 

Taking only the first term of the recursion relation (12.331) . we can approximate the concate- 
nated block failure rate as 

e'^'^c^e'^'. (2.34) 

In the case that e'^' < e, then (12.341) indicates that the error rate of a concatenated code will 
drop doubly exponentially with L, while the size of the code block grown exponentially. This is, 
in essence, what is called the threshold theorem |ITl[T9l|23|, which holds under various different 
conditions, but for our purposes it suffices to say that there is an limitless supply of fresh qubits, 
that the base erasure rate is independent of the size of the circuit, and that the probability of 
erasures on the different gates are independent. The value of e which gives limi_-^oo e*'^' — 
is called the erasure threshold, and obtaining such a value is the main focus of this thesis. In 
practice, this can be taken as the value of e which gives e'^' < e, as long as we assume that gate 
failures are independent, that the gate failure rate does not depend on the computation size, and 
that fresh qubits can be produced on demand. A much more detailed description of how to obtain 
the recursion relation (12.331) and how to extract the threshold will be given in Chapter |5l 



Chapter 3 



Efficient Linear Optics Quantum 
Computation 



A very brief overview of the efficient linear optics quantum computing proposal by Knill, Laflamme 
and Milburn ifTSll is given, focusing on the behavior instead of implementation details. From that 
description, an error model is derived. 



One of the earliest proposals for a quantum computer, put forth by Chuang and Yamamoto H, 
described a system where qubits were encoded in two photon modes - the so called dual rail 
encoding. In order to differentiate quantum states representing photon number states and quan- 
tum states representing the qubit, we follow the convention of using the standard font for number 
states, and the gothic font for qubits, so that the dual rail encoding of the ith qubit would be 
represented by 



where Qi and are the modes corresponding to qubit 1. The main motivation of using photons 
to encode qubits is that photons do not decohere as easily as most other physical systems used to 



3.1 The Knill-Laflamme-Milburn proposal 




(3.1b) 



(3.1a) 
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implement qubits, simply because photons can easily be made to not interact strongly with the 
environment. 

One can use very simple optical elements to perform single qubit operations, and this was 
constructively proven well before the Chuang-Yamamoto proposal ll24l . The elements used are 
called passive linear optics elements, and they are comprised of beam splitters (partially reflective 
mirrors) and phase shifters (delays). These elements have the property that they preserve the 
number of photons, and their behavior can easily be described by how they transform the photon 
creation operators of the modes involved. 

However, it is very hard to make a universal set of quantum gates. In any universal set 
there is at least one entangling gate that requires the interaction of different qubits, and if the 
qubits are encoded as photons, it is very hard to construct such gates without using non-linear 
media to mediate the interaction. In H, a Kerr non-linear medium was proposed to construct 
entangling gates, but Kerr media are notorious for having very high loss. Even measurements 
that require an implicit interaction of the qubits, such as Bell-basis measurements, are impossible 
to perform without failure ll2l1l . It was thought that these facts comprised an informal "no-go" 
theorem for linear optics quantum computing, and various proposals that attempted to build 
quantum computers out of linear optics were shown to require an exponential amount of physical 
resources. 

When Gottesman and Chuang ifTTl demonstrated that fault-tolerant universal quantum com- 
putation could be performed by using quantum teleportation and state preparation, this picture 
changed. The state preparation required for this gate construction scheme depends only on the 
gate being implemented, not on the inputs to the gate, so that if there is a way to prepare these 
states offline, linear optics quantum computation can indeed be performed (a fact that was pointed 
out in dl). One just needs to keep trying to prepare this states until the preparation succeeds, 
or maintain many copies of the successfully prepared states. 

Building on these ideas, Knill, Laflamme and Milburn put forth what is now commonly 
referred to as the KLM proposal for efficient linear optics quantum computation ifTSl . This 
proposal uses single photon sources and detectors, linear optics, and post-selection based on 
measurement outcomes. 

There are two crucial parts of this proposal: the state preparation, and efficient teleportation. 
It was shown that one can build a non-deterministic non-linear sign change operation on a single 
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photon mode, with knowledge of success or failure, using only linear optics and measurements. 
This gate, called NS_i, in effect performs the transformation 

10) + 11) + 12) ^10) + 11) -12), (3.2) 

with a finite probability of success, which is detected by measuring some ancillary modes in the 
gate construction. Loose bounds have been placed on the success probability of constructions 
of this gate using only linear optics 111511 . and it is known that these gates cannot succeed with a 
probability of j or more. One can construct a controlled phase flip, also known as a controlled 
sign or CSIGN, using applications of the NS_i and linear optics elements on the dual rail en- 
coding. Direct computation with these gates is not scalable, but using the teleportation scheme 
proposed by ifTTl . scalability is achieved by using these gates only in state preparation, where 
numerous attempts can be made until the gate performs the desired operation. 

There is still the problem of how teleportation is performed. In the standard teleportation 
protocol, one needs entangled qubit pairs, which can be prepared offline, as well as Bell mea- 
surements. We have mentioned already that ||2T1 showed that such measurements are impossible 
to perform without probabilistic failure - in fact, one cannot distinguish between two of the 
states in the four state Bell basis, and because of that, the probability of failure cannot be made 
less than 50% if the input states are chosen uniformly over the Hilbert space. The alternative 
proposed in ifTSl is to use a modified protocol that relies on the preparation of a larger ancillary 
entangled state, the application of a Fourier transform involving the data to be teleported and 
the prepared state, and measurement of the ancillary states - the Fourier transform can be im- 
plemented efficiently using only linear optics. Given an ancillary state consisting of N qubits, 
this teleportation protocol succeeds with probability a great improvement over the standard 
teleportation protocol. Note that N = 1 is the base case with j probability of success. This 
means one qubit is used for the teleportation, which seems to disagree with the knowledge that 
a Bell state is necessary for teleportation. However, only two of the four modes needed for the 
two qubit CSIGN operation interact with the NS_i gates, so those are the only modes that need 
to be teleported, and each one requires two modes in a Bell state, resulting in one qubit per mode 
teleported. 



Efficient Linear Optics Quantum Computation 



25 



3.2 The Error Model 

Errors are introduced into computation in efficient linear optics quantum computation through 
various sources: failure of CSIGN, photon loss, and finite accuracy of phase- shifters and beam- 
splitters, etc. Here we consider only failures during the application of a CSIGN that can be 
detected by measurement of the ancillary modes. In effect, this is not a detection of the failures 
of the NS 1 gates, but instead a detection of failures that occur during the teleportation of the 
qubits that realizes the CSIGN given the appropriate ancillary state. The NS_i gates are used 
only in the state preparation for the teleportation, so we can choose to simply not use states 
resulting from failures of these gates. 

Two failure modes are discussed here: teleportation failures that are inherent in the protocol 
because of limitations in linear optics, and failures due to photon loss at the detectors during 
teleportation. 

3.2.1 Error model for ideal hardware 

Assuming that all linear optics elements, all detectors and all sources are perfect, failures are 
still possible in the JCLM proposal. This is due to the fact that teleportation succeeds with a 
probability less than unity, although these failures are always detected. 

Consider how the CSIGN is teleported, abstracting from the details of linear optics, as de- 
picted in Figure im The basic idea is to consider the teleportation of two qubits, followed by the 
application of a C S I GN between them. The commutation relations of C S I GN and the gates used 
for teleportation are used to rewrite the circuit in such a way that the CS IGN is applied between 
the prepared states used for the teleportation - see Section l2.4.2l 

In general, the corrections dependent on the measurement outcomes must be applied to both 
qubits being teleported. However, it has been shown that the state preparation can be modified so 
that, if the teleportation of either qubit fails, it affects only that qubit. The protocol can be further 
engineered in a manner such that the failure can be taken to have occurred after the teleportation 
ifTHI . The effect of this type of failure, which is always detected, is that of a Z measurement 
of the qubit that is being teleported, with the measurement outcome made evident through the 
measurement part of the teleportation protocol. This occurs with probability -j;^, where N is 
the size of the ancilla state, which in the case of N = 2 is simply the Bell state (|oo) + |ii)). 
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|a> - 
offline 
Boo) 

|Boo) 




(a) Teleportation followed by a CSIGN. 



|a> 




IP> 



(b) CSIGN commuted backward into the pre- 
pared state, which is equivalent to (a). 



Figure 3.1: The teleportation of a CSIGN gate ifTTlfTHI . The states |Boo) are, in the simplest 
case, the Bell pair -J^ (|oo) + The two qubit measurement is performed over the Bell 

basis. Different Pauli corrections are applied to the output bit depending on the measurement 
outcomes, and this is indicated by the wider lines connected to the Pauli gates. 

One could, in principle, simply ignore the measurement outcome and take the teleportation 
failure to be a Z erasure, because the distinction between the two is classical information and 
classical processing, which we take to be perfect. Formally, this can be shown by the following 
proposition. 

Proposition 1. A projective Z measurement of unknown outcome at a known location is equiva- 
lent to a phase erasure. 

Proof. Expand the probabilistic application of the two possible projection operators for the 
eigenvectors of Z 

Z+ = ^(I + Z) (3.3a) 

Z = ^(I-Z), (3.3b) 

onto a density matrix p. The probability of the outcome being the +1 eigenstate of Z is denoted 
by Pr(Z+), and for the -1 eigenstate is Pr(Z_), and since Z+ + Z_ = I, Pr(Z+) + Pr(Z_) = 1 . 
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Ignoring the outcome of the measurement, p will be transformed into Z(p) described by 



Z(p) = Z+pZ+ + Z_pZ_ (3.4a) 

= l(I + Z)p(I + Z) + I(I-Z)p(I-Z] (3.4b) 

^(p + ZpZ), (3.4c) 

which is a Z erasure, as claimed. □ 



In this case, one would need to take the approach outlined in Section !^. 6. 11 and in the case of 
a CSS code, one would need to measure an X stabilizer in order to determine the syndrome of this 
erasure. This requires the application of CNOTs, and since they are constructed from CSIGNs 
conjugated by Hadamards (see Section l2.4.2l) . it is clear that the type of erasure introduced would 
be an X erasure, not a Z erasure. This is of particular importance for codes that can only correct 
Z erasures, like the codes used to obtain the 0.5 threshold ifT^lfTTI . because it would lead to an 
immediate uncorrectable failure. In our case it desirable to avoid introducing different types of 
erasures, in order to simplify the analysis in Chapter|5l 

In order not to introduce another type of erasure unnecessarily, we cannot use the standard 
fault-tolerant stabilizer measurement from the previous chapter. Instead, we extend a technique 
used for correcting Z measurements, which is described in detail in the next chapter. 

For the purposes of describing the error model during computation, it suffices to say that each 
CSIGN has a probability et of performing an unintentional Z measurement at either the control 
or the target qubit, independently ^ For simulation purposes, we can model this as a Z erasure at 
either target or control qubits independently, since from a Z erasure, with perfect detectors, one 
can easily obtain the measurement value without incurring any cost. Thus, with the first qubit 
being the control and the second being the target, the error model will be taken as 

control failure : Z ® I (3.5a) 
target failure : I ® Z, (3.5b) 

were the failures are taken to occur independently with probability ei. Details of how the mea- 
surement outcome can be used to aid the correction step will be given in Section 14.1.31 
'The subscript i standing for "ideal". 
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3.2.2 Error model for lossy detectors 

Still assuming infinite precision in the parameters of the phase- shifters and beam-splitters, we 
can consider the possibility of photon loss due to detector inefficiencies. The dual-rail encoding 
of qubits ensures that as long as a qubit is properly encoded, there is in total a single photon 
between the two modes corresponding to the qubit, and since linear optics preserves the total 
photon number, measurement of ancillae allows for the detection of leakage from the dual-rail 
encoding. 

In the original proposal for efficient linear optics quantum computation, a robust teleportation 
protocol called RTi is described. This protocol has the property that it is able to detect the 
usual teleportation failures as well as photon losses both in the qubits used for teleportation 
(ancilla or data) as well as in the detectors. Like any possible linear optics implementation 
of teleportation, this teleportation circuit has only a finite probability of success even in the 
ideal case, and the ancilla measurement outcomes will indicate this type of failure that is not a 
consequence of photon loss. In the case of photon loss, however, the measurement outcomes will 
be different from the case of ideal failure, and as a side effect the modes of the qubit at the output 
of the teleportation are replaced with a fresh dual-rail encoded qubit in a fixed state ifTHI . This 
replacement by a fresh qubit corresponds to a total loss of information about the state of the qubit 
over which teleportation failed. We know that in the case of failure without photon loss we can 
model the output by a phase erasure, as demonstrated previously, but in the case of failure due to 
photon loss, we have the ingredients of a full erasure: knowledge of where the failure occurred 
(through the outcome of the ancilla measurements), and complete loss of knowledge about the 
state of the qubit. 

We have already seen that in the case of a phase erasure, the qubit teleportations realizing 
the CSIGN are affected independently - if the teleportation of the control qubit fails, it does 
not imply that the target qubit teleportation fail, and vice- versa. This is not the case for the full 
erasure type of failure. Consider the teleportation of two quantum states, and the subsequent 
application of a CSIGN as in Figure EU Taking a worst case approach, photon loss in the 
control part of the teleported CSIGN (top half of the figure) induces a Z error at the bottom 
whenever there is an X error at the top as well (this includes Y errors, since it is a combination 
of Z and X errors). A similar effect is observed in a photon loss at target part of the teleported 
CSIGN. Clearly there is a classical correlation between the types of errors at the top qubit and 
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the types of error at the bottom qubit - however, photon loss occurs independently at the top and 
bottom part of the teleported CSIGN. Assuming perfect and instantaneous classical processing 
and communication, we could exploit this correlation in order to reduce the number of syndrome 
measurements that need to be performed. Ideally, because CS IGNs are only applied transversely 
between different encoded blocks, one should choose to perform syndrome measurements on the 
block with fewer erasures, since one is less likely to induce a failure that way. In our worst case 
approximation, however, we ignore this classical correlation between the two encoded blocks of 
data, and we will restrict ourselves to the partial description of the two qubit output. Thus, taking 
the combined output to the CS IGN to be p, and the control and target states to be the result of a 
partial trace over p, that is Pcontroi = tttargetP and Ptarget = tr control P, our error model can be seen 
as each type of failure inducing independent superoperators on the qubits as follows 



control failure 



target failure : 



Pcontrol ( Pcontrol ] 

Ptarget ( Ptarget ) • 

Pcontrol Pcontrol] 

Ptarget ( Ptarget ] • 



(3.6a) 
(3.6b) 



Again, we assume that the probability that a photon loss failure has occurred in either telepor- 
tation is independent of the probability that a failure has occurred in the other teleportation. 
Moreover, we assume that only failures due to photon loss occur in this model, and we assume 
the probability of photon loss in a teleportation is given by under this model^. 



3.2.3 A mixed model 

In general, depending on the gate construction, the probabilities of ideal teleportation failure 
and photon loss determines the error operator for failure of the qubit. However, in order to 
place a bound on the accuracy threshold for any gate construction, we consider only the end 
cases, Pr(ideal failure] failure] — 1 and Pr(photon loss failure! failure] = 1 , which are the cases 
described above. It is important to emphasize that the two error models are very different - in 
one case we have Z measurements with known outcomes occuring independently at either qubit 
of the C S I GN, while in the other we have, under our worst case approximation, both full and Z 
^The subscript I standing for "lossy." 
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erasure occuring at either qubits of the CSIGN. The actual error model will be a probabilistic 
mixture of the two models considered here. 

We can restrict ourselves to considering only the end cases because, as will be demonstrated 
in later sections, the erasures that follow from photon loss have a higher cost than the Z mea- 
surements of the ideal model (i.e. more teleported gates need to be applied), so the threshold 
for any mixed error model will fall between the thresholds of these two end cases. This also 
provides further motivation for using essentially the same correction procedure to correct both Z 
measurements and Z erasures, since it simplifies and reduces the circuitry significantly. 



Chapter 4 
Candidate Codes 



In this chapter a comparison is drawn between two small CSS codes, one being the Steane code, 
a [[7, 1 , 3]] code often used in threshold calculations for general errors, and the Grassl code, a 
[[4, 2, 2]] code that is the smallest single erasure correcting code. Brief descriptions of universal 
sets of gates, as well as details of the erasure correction procedure, are given in order to justify 
the preference for the Steane code. 



4.1 Steane Code 



The Steane code Ol |271 was one of the first CSS codes to be discovered. Is it based on the 
[7,4,3] classical Hamming code and its dual, the [7,3,4] Simplex code, yielding a [[7, 1,3]] 
quantum code that allows for very simple fault-tolerant computation. The 7-bit Hamming code 
has parity check matrix 

[11110 

H[7,4,3]= 110 110 , (4.1) 
10 10 10 1 
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so following the CSS construction of self-orthogonal classical codes in Section E.4.11 we find 
that the stabilizer S is generated by 



Mi=X(g)X(g)X(g)X(g) I (g) I (g)I 

M2 = X(8)X(8) I (g) I (g) X (g) X (g) I 
M3=X(g)I(g)X(g)I(g)X(g)I(g)X 

M4 = Z(8)Z®Z®Z®I®I®I 

M5 = Z(g)Z®I®I®Z®Z®I 

M6 = Z(8)I(g)Z®I®Z®I®Z, 



Because the structure of the stabilizer generators is directly derived from the parity check matrix, 
it immediately follows that the derived CSS code has the same minimum distance as the classical 
Hamming code, as argued in Section lT.4.11 

4.1.1 Fault-tolerant Universal Gates 

Being a self-orthogonal CSS code, the Steane code has a very simple fault-tolerant implementa- 
tion of encoded gates. Shor fIE\ was the first to demonstrate how universal fault-tolerant compu- 
tation could be performed on the Steane code by explicitly constructing a fault-tolerant encoded 
Toffoli using only measurements and encoded Clifford gate operations. We follow his approach 
by first giving the encoded Clifford gates demonstrated by Gottesman ifTOl . and then use in- 
sights by Zhou et al. ll33l to demonstrate how a generating set for the Clifford group Gz can be 
constructed. 

First, we need to implement the Pauli operators. Following the stabilizer formalism, we can 
choose members of the normalizer N (S) of the stabilizer that obey the commutation relations for 
the Pauli operators, namely: 



This choice of encoded Pauli operators is equivalent to choosing the encoded basis states, and it 
is equally non-unique. One such choice is: 



{X,Z} = 



(4.3a) 
(4.3b) 



XZ = Y. 



X 



I ® I ® I ® I ® X ® X (g) X 



(4.4a) 
(4.4b) 



Z 



I(g)I(g)I(g)I(g)Z(g)Z(g)Z, 
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but one could just as well call the first operator Z and the second one X, and the computational 
basis would be changed to whatever the eigenbasis of Z is. These operations do not require 
interaction between different qubits so they are automatically fault-tolerant, and this is always 
the case for encoded Pauli operations based on the normalizer of the stabilizer code. 

The next step is to show the construction of the encoded Clifford gates not in V]c- We have 
seen that the Hadamard gate maps between X and Z and leaves the other Pauli operators invariant, 
and, by construction of the Steane code, every stabilizer operator that is made up of only Xs has 
a counterpart that is made up of only Zs, so applying the the Hadamard trans versally preserves 
the stabilizer, and therefore it preserves the codespace. Given our choice of X and Z, it is also 
clear that this operation is the encoded Hadamard itself. The phase gate P can similarly be 
applied qubitwise to map the all Xs stabilizer generators into all Ys operators with an added 
phase factor, but this is only a valid encoded operation if these operators made up of only Ys are 
in the stabilizer. Given that 

PXP^ = lY (4.5a) 
PZPt = Z, (4.5b) 

we would like the weight of the all Xs and all Zs stabilizer generators be a multiple of four so that 
the J complex phases add up to a phase. Fortunately, the Steane code is a doubly-even CSS 
code, so all stabilizer generators have weight four. The operation that this qubitwise P performs 
on the encoded operators is 

P®^(P^)®^ = I®I®I®I®iY®iY®iY (4.6a) 

= -iI(g)I(g)I(g)I(g)Y(g)Y(g)Y (4.6b) 

= -iY = -iXZ (4.6c) 

= P^XP, (4.6d) 

The qubitwise phase gate realizes the encoded inverse phase gate. In order to get the encoded 
phase gate itself, one only needs to apply the inverse phase gate qubitwise. 

These constructions of the single qubit Clifford group operations for the Steane code are 
much more robust and desirable than the constructions for the same gates for the Grassl code 
derivatives, as will be shown later. The main reason for this robustness is the fact that no two 
qubit interactions were necessary for any of these constructions, only single qubit operations. 
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(a) Hadamard Gate 



(b) phase gate 



Figure 4. 1 : Transversal one qubit Clifford gates for the Steane code. 

In the context of the KLM proposal for linear optics quantum computing, this is a powerful 
advantage, because it allows us to take these operations to be erasure free in the ideal model, and 
loss free in our simple lossy model. 

All we need to complete the generating gates of 62 is the encoded CNOT. For the construction 
of the encoded CNOT itself, it is unavoidable that two qubit interactions be present. All CSS 
codes allow for transversal CNOT application as an encoded operation, and in all these cases the 
encoded operation that is applied is itself the encoded CNOT ifTOl . This follows again from the 
fact that the stabilizer generators are made up of either Xs or Zs but not both, and following how 
the CNOT maps between Pauli operators, it becomes clear that the mapping will preserve S x 

'We need to consider mappings that preserve the set S x S because the CNOT is a two qubit mapping, while S 
describes the encoding of only one qubit. 
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Figure 4.2: Encoded CNOT for the Steane code. 
A quick calculation demonstrates that 



Z(g)I 



I®X 



(4.7a) 
(4.7b) 
(4.7c) 
(4.7d) 



as claimed. 

The Clifford group is not sufficient for universal quantum computation. However, almost 
any additional gate that lies outside Qz gives us a universal set, and the gates that are usually 
considered are the Toffoli gate^ or the ^ gate^, which are known to give a universal set when 
combined with the Clifford group. The first fault-tolerant construction of the Toffoli for the 
Steane code was given by Shor ll26l . and a fault-tolerant construction of the ^ gate was given by 
Boykin et al. ll22ll . Neither of these constructions was seen to follow the same general framework, 
but Zhou et al. ll33ll demonstrated how they follow directly from an extension of the work by 



^This is the gate mapping |a) |b) |c) — > |q) |b) |c © ab). It is also known as the controlled-controlled-not. 
^This is the gate that leaves the |o) unchanged and maps |i) e^T It is called a ? gate for historical reasons. 
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H 



PX 



-a 



— Iy> 



(a) 2 gate 



(b) Toffoli gate 



Figure 4.3: Fault-tolerant gates needed for universal computation, as depicted in ll33l . All mea- 
surements are in the Z eigenbasis. The dashed boxes represent states that can be prepared off-line 
fault-tolerantly, the double lines represent classical control dependent on the measurement out- 



come. 



Gottesman and Chuang ifTTl on universal quantum computation using teleportation and single 
qubit operations. 

The detailed analysis of how these constructions were obtained will not be given here, but 
for illustration the resulting circuits implementing the ^ gate and the Toffoli gate are given in 
Figure |431 Zhou et al. show how to prepare the state in the dotted box fault-tolerantly The 
construction of the Toffoli involves not only the three required CNOTs used to perform a sim- 
plified teleportation protocol, but also classically controlled application of CNOTs and CSIGNs 
depending on the measurement outcomes. This Toffoli construction still has a threshold very 
close to the Clifford gate threshold when the state preparation yield similar failure probabilities 
as the Clifford gate failure probabilities 

Given our error model, where CNOTs/CSIGNs are costly compared to single qubit gates, 
it can be intuitively seen that the encoded ^ gate should give a threshold much closer to the 
Clifford gate threshold under the same conditions because it requires fewer applications of such 
gates. This will be the operating assumption for the rest of the analysis, and from now on we 

^Fault-tolerant state preparation under an erasure model is also easier than in a general error model, since it is 
always clear when an error might have occurred 
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focus only on the threshold for Clifford gates. 

4.1.2 Erasure Correction 

The Steane code is a [[7, 1 , 3]] quantum code, so immediately it is clear that it will be able to 
correct all weight two and one erasure patterns. This can be done by measuring the generators 
of the stabilizer of the code and inferring the erasure operators from the commutation properties 
measured - the collection of all measurements is known as the syndrome of the error. 

One can check explicitly which erasure patterns are correctable by applying the modified 
Knill-Laflamme condition to all possible erasure operators for a given erasure pattern. Taking 
the approach of checking the correctability of X erasure and Z erasures separately simplifies the 
analysis somewhat - we can say that a pattern is correctable if the X erasure patterns and the Z 
erasure pattern are correctable. Through this brute force approach, one finds that only 28 of the 
possible (3) = 35 weight three patterns, | of the total, are correctable. 

All stabilizer generators have weight four, so any stabilizer measurement can lead at most 
to failures in four qubits. If we would like to measure the syndrome of a single erasure, we 
simply choose stabilizers which act non-trivially over the qubit that has been erased. If there are 
two erasures, it turns out that one can always choose a stabilizer that will act non-trivially on 
either of the erasures but not on the other one ll32ll . If we consider single failures on any of the 
erasure free qubits being touched by the stabilizer measurement, we find that | of all reachable 
weight three erasure patterns are correctable, but the other ^ are uncorrectable. This is always 
the case, regardless of which of the two erasures we choose to correct, or which stabilizer we 
choose to measure (given the constraint that we do want to correct one erasure when we measure 
a stabilizer operator). 

Depending on the type of erasure, a different number of stabilizer operators need to be mea- 
sured. In particular, since a Z erasure just needs to be checked for commutation against an X 
stabilizer, and a X erasure just against a Z stabilizer, it follows that a full erasure needs two 
stabilizer measurements. 

Exactly how the correctability of these erasure patterns can be easily determined, as well as 
which patterns are reachable after the erasure correcting procedure is applied to a given pattern, 
will be elucidated in Chapter |5l 
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4.1.3 Z Measurement Correction 

As mentioned in the previous chapter, we would like to avoid using CNOTs to measure the X 
stabilizer of a CSS code. An alternative method has been developed that allows for measurement 
of the syndrome and correction of the error to be done in a single step, as long as the outcome of 
the Z measurement is available ifTilfTTlfTSll . 

In order to correct the Z measurement, we need to, indirectly, measure an X stabilizer that 
acts non-trivially over a single Z measurement. Since the Steane code only has weight four X 
stabilizers, we can consider only the four qubits over which the stabilizer acts non-trivially, so 
the stabilizer can be rewritten 

(4.8) 

lyi) (4.9) 

Wi) (4.10) 

Its) (4.11) 

174). (4.12) 

Without loss of generality, take the first qubit to be the qubit that has undergone a Z measurement, 
and we apply an X on all remaining qubits if the outcome is and do nothing otherwise. If 
initially we have some state 

4 

}^at|y,), (4.13) 
then after the measurement (and possible X applications) we will have the state 

4 

|p)=^at|yi), (4.14) 



X (g) X ® X (g) X, 

with +1 eigenvalue eigenvectors 

-^(loooo) + lllll)) = 
-^(loioi) + jioio)) = 
-^(looii) + IllOO)) = 
-^(loiio) + llOOl)) = 
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offline. 



3: 



(Y®Z' 



90° 



T 
i 



90° 



X 



Figure 4.4: Modified teleportation protocol introducing only Z measurements in case of the 
failure of the (Z ® Y)9o° gate 



where, if we consider only the non-measured qubits, we have 

Iti) = \ooo) (4.15) 

172) = lioi) (4.16) 

Its) = Ion) (4.17) 

174) = liio). (4.18) 

In order to obtain the original state back, we need a fresh qubit in the state |+) = (|o) + |i)], 
so we apply the circuit in Figure 14.51 Since this circuit is made up of only Qz gates, we can 
consider the teleportation of the state |p) before applying Figure l43l and then simply commute 
that part of the circuit into the state preparation part of the teleportation protocol - the only qubit 
introduced is in the fixed state |+), so it can be taken into the state preparation part as well. 

It is important that we use the slightly modified teleportation protocol shown in Figure l4!4l 
In this qubit teleportation protocol, the only source of failure is the two qubit gate (ZY)9oo, 
which is performed through a CSIGN conjugated by one qubit Gj gates (the other two qubit 
gates are performed similarly). This teleportation protocol guarantees that if there is a failure in 
the (ZY)9oo gate, it will translate to a Z measurement of the qubit being teleported. Because of 
the structure of the circuit in Figure l431 only the Z corrections of the teleportation will propagate 
through the CNOTs, and therefore so will the Z measurements. 

If any of the teleportations fail, the corresponding qubit undergoes a Z measurement, but 
also the correction target, the qubit we were trying to recover, cannot be corrected. Taking each 
teleportation to fail with a probability ei, the correction will succeed with probability (1 — ei)^, 
and if we take Z to signify the Z measurement, we have the following transition probabilities on 
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IP) <^ 

Figure 4.5: Circuit taking the state |p) into an eigenvalue +1 eigenvector of X ® X ® X ® X. 



1+) 



Figure 4.6: Full circuit for correction of a single Z measurement. The CNOTs can be commuted 
backwards in time and combined with the offline state preparation for the teleportations. 
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the four qubits in question 

Pr(I (g)I(g)I(g)I|Z(g)I(g)I(g)I) 
Pr(Z ®Z®I®I|Z®I®I®I) 
Pr(Z ®Z®Z®I|Z®I®I®I) 
PrfZ ®Z®Z®Z|Z®I®I®I) 



ei(1 - ei) 



(4.19) 
(4.20) 
(4.21) 
(4.22) 



with the other patterns with equal weights being obtained by permutation, preserving the proba- 
bilities of transition. 

In this analysis, we have neglected the remaining three qubits of the Steane code, and we 
have assumed that the state described by the four qubits in question was a pure state, which is 
not true. In reality, the state of those 4 qubits will be a mixture, because the seven qubits of the 
code block are entangled, but similar results in that case follow by linearity from the results given 
here. 

With the method described here, we can correct Z measurements without directly measuring 
X stabilizers. We can employ the exact same method to correct Z erasures by first measuring 
the erased qubit in the Z basis, and then simply replacing it with a qubit corresponding to the 
measurement outcome. Once the measurement outcome is known, this same procedure can be 
used. 

In the case of a qubit loss in any of the three teleportations, we take the worst case approach 
once more. Thus, the qubit being teleported is completely erased, and the Z erased qubit that 
is being corrected is Z erased once more. Each of the three teleportations can lose a photon 
independently, so the transition probabilities can easily be computed for the lossy model as well, 
yielding 

Ci? (4.23) 
(4.24) 

- eif (4.25) 

- e,) (4.26) 
(4.27) 

with the factors involving 6, the probability of a detector failing, arising due to the fact that we 
need to perform a Z measurement explicitly. 



Pr(I (g) I ® I ® I|Z ® ] 




= (1 


-5)(1- 


Pr(E ® I ® I ® I|Z ® ] 




= 6 




Pr(Z O E ® I ® I|Z ® ] 




= (1 


-6)e,(1 


Pr(Z ® E O E ® I|Z ® ] 




= (1 




Pr(Z ® E ® E ® E|Z ® ] 




= (1 


-5)ef, 
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4.1.4 The Encoded Error Model 

In the case of a failure in correction an error in a code block, it is important to understand what 
kind of error is induced in the encoded qubit. This encoded error model is crucial for predicting 
the performance when the information is encoded multiple times, an approach that will be used 
in the Chapter |5l 

For the error models considered in Chapter|21 this is a relatively simple problem that can be 
solved with the help of the stabilizer and normalizer formahsms. 

Error Model for Ideal Hardware 

In this error model, all erasures are Z measurements, and from the structure of the Steane code, 
it is known which patterns of erasures are correctable and which ones are not. Consider the case 
of the Z measurement of the last 3 qubits of a code block, which is not correctable. Since the 
encoded Z operation consists of Zs applied to the last 3 qubits, as described by (I4.4bl) . it is clear 
that this uncorrectable error is equivalent to an encoded Z measurement on the encoded qubit. 
Consider, on the other hand, the case of Z measurements on all qubits of the code block. Since, 
given the stabilizer operator M4 from (lO) . ZM4 = Z(g)Z(g)Z(g)Z(g)Z(g)Z(g)Z, itis clear 
that this is also equivalent to a Z measurement. It turns out that because of this, all uncorrectable 
weight three measurement patterns (as well as the weight seven pattern) are equivalent to a Z 
measurement. All uncorrectable measurement patterns of intermediate weight can be thought of 
as a uncorrectable weight three measurement pattern followed by some additional measurements, 
so that all of them are equivalent to a Z measurement. 

In practice, it is also important to consider the cases of correctable erasure patterns, since a 
code block may end up in such a state after some finite number of error correcting rounds. This 
is extremely simple for the case of ideal hardware, since we can simply measure any number 
of qubits perfectly, allowing us to force an uncorrectable measurement pattern, which as shown 
above is always a Z measurement. 

Thus, the encoded error model is identical to the base error model up to different probability 
distributions. 
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Error Model for Lossy Detectors 

Given the correspondence between Z erasures and Z measurements of unknown outcome^, fail- 
ures consisting only of Z erasures have a behavior that parallels failures consisting only of Z 
measurements. 

The stabilizer generators for the Steane code that are tensor products of Xs only have exactly 
the same structure as the stabilizer generators that are tensor products of Zs only, and a full 
erasure is the composition of a Z erasure and an X erasure^, so uncorrectable erasure patterns 
consisting only of full erasures are equivalent to encoded full erasures. 

Uncorrectable erasure patterns that include both Z and full erasures need to be considered 
as composition of Z and X erasures patterns. With the error correction scheme described in the 
previous sections of this chapter, we can only have X erasures on qubits that have Z erasures as 
well - that is to say that we only have Z or full erasures, but no X erasures. If the X erasure 
pattern is uncorrectable, then the encoded failure will be an encoded full erasure, since we are 
guaranteed to also have an uncorrectable Z erasure pattern. On the other hand, if the X erasure 
pattern is correctable, we can simply attempt to correct X erasures until either there are no more 
full erasures (in which case we would only have the uncorrectable Z erasure pattern, which 
corresponding to an encoded Z erasure), or until the X erasure pattern is uncorrectable (in which 
case we would have the encoded full erasure). 

Correctable erasure patterns are a little more complex since they require the measurement 
of all qubits in order to be interpreted as encoded errors, and because of the lossy nature of the 
detectors, the outcome could be either an encoded Z measurement, or an encoded full erasure. 
The probabilities of each of these outcomes depends on the probability of failure of the qubit 
measurements, so this extra measurement needs to be regarded as a separate step after the desired 
rounds of error correction. In the best case, where no measurements fail, we obtain an encoded 
Z measurement, while in the worst case, we obtain an encoded full erasure. As shown before, 
measuring all qubits in the Z eigenbasis is equivalent to an encoded Z measurement. Moreover, 
because of the structure of the Steane code, even though some of these measurements may fail, 
the encoded measurement outcome may still be inferred by applying classical error correction 
techniques to the outcomes. If recovery is not possible because of the high number of errors, 

5 See (|333 in Chapter|3l 

^The order of which partial erasure occurs first is irrelevant. 
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then the failure is an encoded full erasure. We can take these correctable erasure patterns to be 
encoded Z erasures, since, like Z erasures, they require Z measurements to be performed as the 
first step of the error correction, and failure of the measurement yields a full erasure. 

Once again, the encoded error model is identical to the base error model up to a different 
probability distribution. 

A Mixed Error Model 

In a mixed error model, where we have both the limitations of linear optics and the post-selection 
based construction of the CSIGN as well as the imperfect photon detectors, the error model 
changes only quantitatively. At higher levels of encoding we still expect to find Z measurements 
as well as Z and full erasures, but qualitative analysis shows that there is a bias towards the more 
benign error model. 

A failure pattern than includes both Z erasures and Z measurements (but no full erasures) 
can be taken as an encoded Z erasure at no additional cost, or may be taken as an encoded Z 
measurement at the cost of measuring the Z erasures, possibly introducing full erasures which 
may lead to encoded full erasures. If the measurement can be inferred without measuring the Z 
erasures by using classical error correction, the probability of encoded full erasures is reduced. 
Failure patterns that include full erasures are treated similarly. 

The procedure for turning correctable failures into failures at a higher level of encoding could 
either yield Z measurements, full erasures or Z erasures. In particular, there is no need to explic- 
itly measure qubits that have been affected by an unintentional measurement, which also reduces 
the probability of encoded full erasures. 

Both these opportunities for gain over the purely lossy error model indicate that the mixed 
error model should always allow for improvement over the purely lossy error model, and that it 
will be naturally biased towards the ideal detector error model. 

Worst Case Analysis 

For the calculations made in Chapters |5l we take all encoded errors in the error model for lossy 
detectors to be full erasures. This worst case approach simplifies the calculations significantly, 
and although it gives a looser lower bound to the error threshold, it will be shown to be enough 
to match predictions in ifTTlfTHI . 
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Formally one usually takes the threshold to be the break even point between the encoded and 
the base error rates. Here we will take the threshold to be the point where the encoded error 
rate is half of the base error rate. This is because the probability of full erasure of a qubit given 
a failure has occured is the same as the probability of a Z erasure given a failure has occured 
- that is, both occur with conditional probability ^. Thus, since we take all encoded failures 
to be full erasures, we apply the break even condition to full erasures only, which translates to 
the requirement that the overall encoded failure rate be half of the base failure rate. When the 
threshold is explicitly calculated in Section |5?T1 this relationship will be formalized more clearly. 



4.2 Grassl Code 

Grassl, Beth and Pellizzari lfT2l have shown that the shortest possible erasure code is a [[4, 2, 2]] 
code, which is referred to here as the Grassl code. This is a CSS code derived from the classical 
code consisting of all even weight strings of length four. It corrects a single erasure, the general 
class of such distance two codes having been discussed in detail in flOl . The stabilizer group S 
of this code is generated by the operators 

Mi=X®X®X®X 

(4.28a) 

M2 = Z®Z(g)Z(g)Z, 
and a generating set for the encoded Pauli operators is 

Xi=X(g)X(g)I(g) I 

X2= X ® I ® I ® X ^^29a) 
Zi=Z(g) I (g)I(g)Z 
Z2=Z(8)Z(g)I®I, 

The fact that all operators in the generating set of the encoded Pauli operators are of weight 
two proves that this is a distance two code - that is, there are weight two error operators that 
will map valid codewords into different valid codewords, and will therefore be undetectable and 
uncorrectable. 
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4.2.1 Fault-tolerant Universal Gates 

One difference between the Grassl code and the Steane code that is immediately obvious is the 
fact that the Grassl code encodes two qubits. This raises the question of how to concatenate 
the code, and how the different choices affect the performance. Even if this issue is resolved, 
because we have two encoded qubits, we must consider different types of encoded operations 
dealing with the different encoded qubits. While there are very simple constructions of some of 
these gates - such as an encoded CNOT between the two encoded qubits in a block or the SWAP 
between these same qubits - the construction of other types of CNOTs have a very high cost 
because of the use of multiple rounds of teleportation and multiple ancilla preparations - some 
gates that suffer this problem are the encoded CNOT between the individual encoded qubits of 
two different blocks. 

On the other hand, one may choose to bypass this problem by modifying the Grassl code to 
encode a single qubit. That can be done by choosing one of the encoded Pauli operators over 
a single encoded qubit and making it a stabilizer. Any such choice is valid since we know that 
all operators in N(S) commute with the elements of S (see Section 12.4.21) . and that results in 
a stabilizer code with three stabilizer generators, as desired (n — k = 4 — 1 since we want a 
[[4, 1]] code). It is easy to check that the minimum distance of the code remains the same, so we 
are able to correct a single erasure with this code. The main problem is that the resulting code 
is not a self-orthogonal CSS code, but a more general form of CSS code. While the encoded 
CNOT can easily be implemented by qubitwise CNOTs, neither the encoded Hadamard nor the 
encoded Phase gates can be implemented as easily. Instead, it is necessary to use ancillary states, 
multiple CNOTs and single qubit operations to generate these two remaining encoded Clifford 
gates ifTOl . This is highly undesirable for the error models under consideration because this 
translates to encoded single qubit operations that can introduce erasures, while in the Steane 
code those operations can never introduce erasures. 

A similar problem is found even if the Grassl code is maintained in its original form. The 
encoded single qubit operations that are applied only to one encoded qubit, while the identity is 
applied to the other, require complex protocols involving teleportations and measurements. Thus, 
there is strong evidence that the Grassl code is suboptimal for fault-tolerant quantum computation 
in the error model given for linear optics quantum computation. 



Chapter 5 

The Erasure Threshold 



A novel method for determining the recursion relation for the encoded failure rate is given. This 
method describes the erasure correction procedure as a random walk over equivalence classes 
of erasure patterns, allowing for a compact description that can give the recursion relation for 
any number of erasure correction attempts. Using this technique, the recursion relations for the 
Steane code under the ideal and lossy error models are given, and values for the threshold are 
estimated numerically. 

5.1 Calculation of the Threshold 

In order to calculate the error threshold, we need to calculate the recursion relation for the proba- 
bility of error e''^' in encoding level L in terms of the probability of error e''^^^' in encoding level 
L — 1 . For simplicity, we consider only the relationship between the first level of encoding and 
the base error rate, that is, 

oo 

e(^' = Y_ Cne^ (5.1) 

rL=Ta 

where e is the error rate at the physical qubits, Cn are constant integers dependent only on n, and 
ra is the minimum number of failures causing an unrecoverable error, as we have seen in Chapter 
11 

In the general case one needs to keep track of the possible propagation of undetected errors, 
but the erasure error model allows us to ignore this complication because the presence and lo- 
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cation of erasures is known by definition. Thus, we only apply the erasure correction procedure 
when we know an erasure has occurred, and any additional failure that occurs during erasure cor- 
rection is immediately flagged, so we can keep applying the erasure correction procedure until 
the data is erasure free or we find the resulting threshold for the number of trials is acceptable. 

Since the Steane code is a [[7, 1 , 3]] code, we know it can correct any one or two erasure 
patterns on any given code block, so it is impossible that only one or two erasures at one level 
will cause an erasure at a higher level of encoding. Some weight three erasure patterns are not 
correctable, so it is possible for three erasures at one level to cause a failure at a higher level. 
Therefore, the leading term of the series expansion of (15.11) will be cubic power in e, i.e. 



We cannot restrict the analysis to erasures that occur during computational steps only, but we 
also need to consider erasures that are introduced during the erasure correction steps, since both 
are performed in the same physical system. 

In essence, we need to consider erasures that occur during 

• computation 

• interaction with ancillae used for syndrome measurement 

• measurement of ancillae used for syndrome measurement 

Once the recursion relation has been computed, the threshold is found by solving an equation 
describing the break even point, that is, the probability e at which the encoded failure rate is 
equal to the base failure rate 



whichunder our simplifying assumption, gives, for e' ' < e the required condition ItmL^ooe = 
0. Note that for the lossy detector error model, we modify the break even condition to be the worst 
case break even condition, since we take all failures to lead to encoded full erasures. Moreover, 
the probability of a single qubit full erasure is half of the total probability of failure, yielding the 
worst case break even condition 



oo 




(5.2) 



n=3 




(5.3) 




(5.4) 
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Because the recursion relation describing ' is of very large (if not infinite) order in e, these 
equations cannot be solved exactly, but it can be approximated numerically. It is not clear how 
many leading terms of the recursion relation are sufficient to accurately determine the error 
threshold, so we would like to be able to easily compute as many of the terms as possible. 

One approach is to consider all possible correctable erasure patterns. After an application of 
the encoded C S I GN, the only source of failures in our error model, the erasure pattern determined 
by the failures of the individual gate is known. By attempting to correct the erasure, one applies 
CSIGNs between the data and the ancilla in order to measure stabilizer operators (see Section 
12.6.11) . Note that we can restrict ourselves to CSIGNs because a CNOT is simply given by 

CN0Ti,2 = (I ® H)CSIGNi,2(I ® H), (5.5) 

and if we only apply the second Hadamard gate after all erasures have been corrected, we have 
the advantage of having an error model that still consists only of E (full erasures) and Z (Z 
erasures). 

Any failure during the erasure correction will change the erasure pattern in the data into some 
other erasure pattern. Since we know the probability of the different failures occurring, and we 
know which operators will have to be measured in order to correct a given erasure pattern, we 
may describe the process by transition probabilities between different erasure patterns. In the 
case that the erasure rate at each time step is independent of the previous time steps, what ends up 
being described is a Markov chain, or a random walk in a graph. Therefore, given the set of pos- 
sible correctable erasure patterns Wi, we want to calculate all transition probabilities Pr(Wi|Wj) 
for some error model. If we have a description of the initial distribution of erasure patterns, we 
can compute what is the probability of each erasure pattern after an erasure correcting step by 
following a simple application of Bayes' theorem 

Pr(Wt,tt+i) = }^Pr(Wt|Wj)Pr(Wj,tO, (5.6) 

where we take tt+i to be the time instant immediately following tt. Once the transition proba- 
bilities have been computed, by repeated application of this formula one obtains the probability 
distribution of all erasure patterns at time In, which corresponds to the distribution after N 
rounds of the erasure correction procedure. At that stage, any erasure pattern that has weight 
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greater than is considered a failure, and if we compute 

£<!)= Y_ Pr(W|,tN), (5.7) 

wt(Wt)>0 

we have the recursion relation for N rounds of erasure correction' . 

We may rewrite the above transition probabilities into a transition matrix !P with elements 

[y]y = Pr(Wt|Wj]. (5.8) 

Then the probability distribution of the erasure patterns may be described by a column vector 
J (to) given by 

P(to)]^=Pr(Wi,to). (5.9) 
This allows for a recursive definition of the probability distribution at time 

a(tN)=ya(tN-l)=rt(to). (5.10) 

Because of the way in which T was defined, it satisfies the conditions of a stochastic matrix, 
which guarantees that J(tN) as defined in (15.101) is a probability distribution as long as J(to) is 
one as well. 

It is also worth noting that certain erasure patterns in the above description have special 
properties. Consider, for example, the erasure pattern of weight zero. There is no need to apply 
the erasure correcting procedure to this erasure pattern, so the probability of transition from this 
pattern to any other of non-zero weight is zero. This is what is called an absorbing state. Other 
examples of absorbing states are uncorrectable erasure patterns. Since we know that there is no 
way we can go from an uncorrectable erasure pattern to a correctable one, there is no point in 
applying the erasure correcting procedure either. 

This approach, although simple to describe, is by no means efficient. The number of possible 
erasure patterns grows exponentially in the number of qubits in the code. 

5.1.1 Equivalence of Error Superoperators 

There is great redundancy in the description of the transition probabilities at each erasure correc- 
tion step. In the context of the Steane code, if we observe all transitions from a weight one erasure 
'Note that as N ^ oo, the probabiHty distribution converges to a stable distribution. 
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Figure 5.1: Portion of a naive Markov chain describing the erasure correction of some weight 
one erasure patterns. Here Pr(succ) indicates the probability of success in correcting one phase 
erasure. 

pattern, we notice they are very similar. Namely, they will have the same number of transitions 
to patterns of any given weight, and each of these transitions will have the same probability. 

The regularity of the graph depicted in Figure ISTTl is no coincidence - it is, in fact, a direct 
consequence of the symmetries of the erasure correction code being employed. Intuitively, many 
of these erasure superoperators are equivalent, and can be 'bundled up' into a single represen- 
tative superoperators, eliminating the redundant edges and vertices in the graph in favor of a 
minimal description that differentiates between the operational behavior of the superoperators - 
that is, their correctability, and their transitions to other superoperators. Formally, we need to de- 
scribe an equivalence relation between erasure superoperators. In order to do that, we introduce 
the concept of a code automorphism. 

Definition 2. The automorphism group Aut( C) of a stabilizer code C is defined as the subgroup 
of the group generated by qubit permutations and qubitwise applications o/H and P leaving the 
stabilizer of C invariant / l25l/ . 

The code automorphism can be used to define the equivalence relation we need. 

Definition 3. For a given quantum code C with automorphism group Aut(C), and given the 
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Pauli group over n qubits Vn, we define the relation =c over Vn such that 



E=cF ^3g e Aut(C) : E = gFg 



(5.11) 



Given this definition, some elementary results about =c can be given. 
Proposition 2. =c is an equivalence relation over Vn- 

Proof. Note that Aut(C) is a group, so, by definition, is it closed under product, it contains the 
identity, and every element has an inverse. In the following, take Et G Vn for any index i. We 
show that =c has each of the properties of equivalence relations: 
Reflexivity: It is always the case that Et =c E|. 

This follows trivially from the fact that, because Aut(C) is a group, I G Aut(C], so g = 
g-i = I and E, = gE,g-\ 

Symmetry: If E^ =c Ej, then Ej =c Et. 

It is given that Et = gEjg^\ so conjugating both sides by g^^ yields g^^E^g = Ej. Since 
g"^ G Aut(C), it follows that Ej =c Ei 

Transitivity: If Et =c Ej and Ej =c Eic, then Et =c E^. 

It is given that Et — gEjg^^ and Et = h,E]^h,^\ By substitution, Et = gh,Ei^h,^^g^\ and 
since g, H G Aut(C), it follows that gh, h^^g^^ = (gH)^^ G Aut(C) and that Et =c E]<. 

This concludes the proof. □ 

Theorem 1. If error operators {Et} satisfy the modified Knill-Lafiamme condition for some code 
C, then so do {gEtg''"}, where g G Aut(C). 

Proof. The first modified Knill-Laflamme condition states that 



for any basis codewords |ct) , |Cra) G C. Since g G Aut(C), so is g ^ — q\ and one can insert 
the identity gg^^ = g^^g = I in the product (ctl Et |ct) 



(ctlEtlct) = (cjEtlcTT^) 



(5.12) 



(ctlEtlct) 



(ctl g"^gEtgTg |ct) 
(ctl gEtg"^ |ct) 



(5.13a) 
(5.13b) 
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where the automorphism maps the basis {|Ci)} into {|Ci)}, and since all code automorphisms are 
unitary, both bases are orthonormal. Moreover, both bases describe the same space - the code 
space - so operators satisfying the Knill-Laflamme conditions in either basis are correctable. 
More explicitly, the first condition can be re- written 

{ci\ gEtg"^ \ci) = (cj gEtg^ |c^) , (5.14) 

so clearly it is satisfied for any {gEig^^} regardless of the automorphism g G Aut(C). 
The second modified Knill-Laflamme condition states that 

(cil Ei |c^) = for {ci\c^) = 0. (5.15) 

Again, it follows quite simply that 

(cilEilCnx) = (cil g^gEig^g |c^^^) (5.16a) 

= (ctl gEtgt |c^) (5.16b) 

= 0, (5.16c) 

and that 

(ci|c^) = (cdg^glcTT^) (5.17a) 

= (cilc^) (5.17b) 

= 0, (5.17c) 

concluding the proof. □ 

Note that, because Aut(C) is made up of permutations and tensor products of single qubit 
unitary operations in N(S), where S is the stabilizer of the code C, the weight of gEtg^ is the 
same as E for any E G Vn and g G Aut(C) for any code C. 

It has been shown in Chapter El that an erasure pattern is a uniform, convex sum over Pauli 
erasure operators. In that case, we can extend the definition of =c to be a relation over erasure 
patterns as well. We say W =c K if and only if there is a g G Aut(C) mapping between the set 
of all erasure operators {E|}| in the convex sum W to the set of all erasure operators {Fjl^ in the 
convex sum K. 



54 



Erasure Thresholds for Linear Optics 



In some cases, however, we want to change which operations g are allowed in this definition 
because the error correction procedure may not be invariant under all code automorphisms. The 
problem becomes clear if we consider an error model that allows for Y erasures or X erasures, 
which under this definition are equivalent to Z erasures. The Y erasure requires two error cor- 
rection steps in the best case and thus has different behaviour from a Z. Similarly, we have used 
different circuitry to correct X erasures and thus, if X erasures were not always accompanied by 
Z erasures, they should not be considered equivalent to Z erasures. For the error model consid- 
ered here, this distinction is not necessary, and the definition in term of the automorphism group 
is suitable and simple to analyze. However, in general, we want the equivalence relation to relate 
erasure patterns that have equivalent error correction behaviour, and a different error model may 
require a different definition for the equivalence relation between the erasure patterns. A simple 
restriction, for example, would be to consider the autopermutation group of the code, a subgroup 
of Aut(C) that restricts g to qubit permutations only. This restriction, however, is not necessary 
for the error models and correction procedures considered here, and for the rest of the thesis, we 
consider =c as a equivalence relation over erasure patterns as defined in the previous paragraph, 
which includes the full automorphism group. 

If, for two correctable erasure patterns W, K we have W =c IK, then failures during the 
erasure correcting procedure of W lead to error operators equivalent (under =c) to the error 
operators resulting in failures during the erasure correcting procedure of K, with the same prob- 
abilities. 

Say we have W, and we measure an operator M G S in order to attempt to correct a sin- 
gle erasure. Each of possible \i failures during the measurement leads to the erasure patterns 
Wi , ■ ■ ■ , with respective probabilities Pr ( Wi |W) , ■ ■ ■ , Pr( W^,| W) . In order to correct K, we 
can apply a code automorphism to map it to W, then measure the same operator M, which leads 
to the same possible erasure patterns under failure, with exactly the same probabilities. Apply- 
ing the inverse of the automorphism leads to equivalent erasure patterns under =c, with the same 
probabilities. Note that the application of these automorphism does not need to correspond to 
operations in the physical system. Permutations simply relabel the qubits, while applications of 
H and P relabel the stabilizers, so that one may simply use this re-labeling to determine which 
operators need to be measured. 

Thus, when analyzing the transition probabilities, we need only consider transitions between 
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equivalence classes under =c, greatly reducing the number of transitions and erasure patterns 
that need to be analyzed individually. 

In practice, these equivalence classes can be found through computer algebra packages such 
as GAP |T| or MAGMA (5\ . In the case of the Steane code, we have already seen some transver- 
sal (qubitwise) unitary operations that preserve the stabilizer (and thus the code space) - namely 
H, A quick search through GAP demonstrates that the qubit permutations^ 



(1,2)(5,6) (5.18a) 

(2.4) (3,5) (5.18b) 
(2,3)(4,6,5,7) (5.18c) 

(4.5) (6,7) (5.18d) 

(4.6) (5,7), (5.18e) 



generate all permutations that preserve the code space. The group generated by all these permu- 
tations and the unitaries described above is the automorphism group of the Steane code. It has 
order 168, and it is isomorphic to the group of invertible 3x3 matrices in GF(2). 

5.1.2 Z Measurement Threshold Calculation 

Focusing our attention on the error model outlined in Section 13.2.11 and on the Markov chain 
method outline in the previous section, we can calculate the recursion relation for the uninten- 
tional Z measurements that can occur during teleportation. 

It has already been shown that there is a code that has a Z measurement threshold of 0.5 for 
this error model ifTTI . but that code can only correct single Z measurements in a block. In Chapter 
131 an alternative method for correcting Z measurements was described, with the advantage that 
only Z measurements are introduced as a result of failure, thus we can consider only the subgroup 
of Aut(C) corresponding to the qubit permutations, which is the code automorphism of the 
classical Simplex code [7, 3,4] used to construct the Steane code. Any two weight one erasure 
patterns are equivalent under the automorphism relation. Similarly, any two weight two erasure 

^In this notation, a permutation (a, b, c) means that qubit a takes the place of b, b takes the place of c, and c 
takes the place of a. Concatenation of two of these permutations is equivalent to applying the rightmost permutation 
first, and applying the permutations, in order, from right to left, just like products of operators. 
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patterns are equivalent as well. This is not the case for weight three erasure patterns, however. 
We know that Z is a weight three Pauli operator consisting of only Zs and identities, so it will 
be an element of the convex sum I®I®IC?)I®Z®ZC?)Z, and thus this erasure pattern will 
be uncorrectable. By applying the permutations of the code automorphisms, we find that only 
a subset of the weight three Pauli operators are generated. That is, the equivalence class which 
includes this pattern does not include all the possible erasure patterns, it only includes 
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(5.19) 



All 28 = (2) — 7 remaining weight three erasure patterns are equivalent to each other. Li order 
to see that they are correctable, consider 

I®I®I®Z®I®Z®Z (5.20) 

We can correct the rightmost erasure by measuring the stabilizer M3 from Section HTTl reducing 
the pattern to a weight two erasure which is always correctable. 

Thus, we are left with four equivalence classes corresponding to correctable erasure patterns 
of weight 0, 1, 2 and 3, and one equivalence class corresponding to uncorrectable erasure pat- 
terns (with weight 3 or more). The transitions between the different equivalence classes can be 
determined by analyzing a single erasure pattern in each class under one erasure correcting step. 
It suffices to say that for each case a weight 4 stabilizer needs to be measured, which requires 
four CNOTs between the data and an ancilla, and the measurement of the ancilla. Failures at any 
of these CNOTs must be considered, which is a simple but tedious process, described in detail in 
Appendix |A1 leading to the simplified pictorial view of Figure ls!2l 

From the Markov chain description one obtains the error rate recursion relation 

' = 56el + 406et + 3878ef - 1 29675ef + 11 6481 5e[ + • ■ • (5.21) 
with a threshold of approximately ei = 0.1 15. 
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Figure 5.2: Graph describing the Z erasure correction procedure for the Steane code under the 
ideal error model. The probabilities associated with each transition are omitted for readability. 



Theoretical Predictions for Ideal Error IVIodel: log-log scale e/ei*'' vs. E| 
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Figure 5.3: Plot of under the ideal error model. If the ratio is greater than 1, the encoded 
measurement probability is lower than the uncoded measurement probability. 
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The correctness of the leading term in ( I5.21I ) can be quickly checked. First of all, because 
we are dealing with a distance 3 code, all measurement patterns of weight up to two can be 
corrected. That immediately sets the exponent of the leading term to three. In order to determine 
the coefficient, we must count how many ways there are to reach an uncorrectable measurement 
pattern with 3 errors during the correction procedure, assuming of course we are doing a minimal 
number of correcting rounds to be able to correct a weight three measurement pattern. First, we 
have seen that ^ = ^ of the weight three patterns are uncorrectable, that leaves the coefficient 
at ^ = 7. Now consider the correction of a weight two pattern: any single additional failure 
brings it back to a weight three pattern, of which j are uncorrectable, so we add (J^ |(^) = 21 . 
In the case of weight one patterns, which are also always correctable, we could have two failures 
at one, or one failure at a time, adding (^) 3 ((2) + (^) (^)) = 28. The total is, as expected, 56. 

5.1.3 Full Erasure Threshold Calculation 

In the lossy error model, things become a little more complicated. In particular, we have a greater 
variety of erasure patterns, which in our case lead to a greater variety of equivalence classes and 
transitions. In order to make it clear, note that the qubitwise operations that are allowed in the 
definition of the automorphism of a code are simply single qubit Clifford operations, which in 
essence permute the labeling between the non-trivial Pauli matrices X, Y, Z. Thus, E is invariant 
under these operations if we consider a basis change given by the single qubit 62 operations, 
but X,Z,Y are not. In fact, they undergo the same label permutations as the Pauli matrices 
themselves. That means that now, instead of keeping track of the weight of the erasure patterns, 
we must keep track of the number of different types of erasures. Some general properties still 
hold, such as being able to correct any weight one or two erasure patterns, and | of all weight 
three patterns, as long as they consist of a single type of partial erasure. This is not as bad as it 
appears at first, we can think of erasure patterns of a single type of partial erasure acting together 
to result in an erasure pattern with different types of erasures. That is, we may consider 



E(g)I®I(g)Z(g)X(g)I(g)I 



(5.22) 
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as a series of erasure patterns 



X(g)I(g)I(g)I(g)X(g)I(g)I 



(5.23a) 



Z (g) I ® I ® Z ® I ® I (g) I 



(5.23b) 



applied to the same code block. So, because any weight two partial erasure is correctable with 
the Steane code, any permutation of (15.221) is correctable, even though it has weight three, so 
in general, we need to keep track of a tuple that has the number of E, X, Z individually, and 
equivalence classes between these tuples can be developed. 

To simplify matters, because we are dealing with a self-orthogonal CSS code, we can restrict 
ourselves with the measurement of Z stabilizers only through the application of CSIGNs, so the 
only types of errors that are introduced according to the lossy model are either E or Z. Following 
a decomposition like (I5.23I) . we find that the Z partial erasure pattern will have a larger weight 
than the X partial erasure pattern if any E occurred in the original pattern. In that case, because 
Zs are just as frequent as Es during computation in our model, it is wiser to correct a Z first. If 
we choose to correct a partial Z without a corresponding partial X, and if we choose to correct 
an X first when there are only E, we ensure that the pattern consists of only E and Z at all times. 
The analysis can then be restricted to 2-tuples describing the number of full erasures and the 
number of Z erasures - graphically, they are described by [m, n] where ra is the number of Es, 
and n is the number of Zs, with the advantage that the weight of the patterns is given by ra + n. 

The same correctability results in the previous section apply here, so we find that all patterns 
of weight one and two are correctable, and that for any equivalence class of patterns with weight 
three, | of the patterns are correctable, and ^ are not. Again, we bundle all uncorrectable patterns 
into a single class. The resulting graph can be seen in Figure l54l and the transition probabilities 
are described in detail in Appendix |B] The main difference between the lossy model and the 
ideal model is that for the lossy model we consider measurement failures, with a probability of 
failure 6, along with a probability of CSIGN failure ei as usual. The transition probabilities are 
given as a function of both considering them to be independent parameters^, but in reality they 
are not. However, this dependence can be recovered by writing one parameter in terms of the 
other in the general expression obtained from the Markov chain description. 

^By independence, I mean algebraic independence, not statistical. 




Figure 5.4: Graph describing the erasure correction procedure for the Steane code under the lossy 
model. Each states is labeled [m, n] , where m is the number of full erasures, and n is the number 
of Z erasures. The probabilities associated with each transition are omitted for readability. 
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For the case where we take 6 = 



' = 350ef + 4739et - 1 2404ef - 355600e^ - 30871 1 0e[ + ■ ■ 




(5.24) 



yielding a threshold of approximately ei, = 0.0324. 
For the case where we take 5 = ei. 



ei^' = 1 050ef + 331 73ef - 46242ef - 6861 701 - 11 8743847e[ + 



1 




(5.25) 



which yields a threshold of approximately et = 0.0]78. This threshold is only valid if the 
encoded qubit measurement failure rate 5'"^^ also vanishes for L ^ oo. Because of the structure of 
the CSS code, where the encoded states are superposition of classical codewords, measurement 
failures can be corrected for. In the case of the Steane code, the classical codewords that are 
measured are codewords in the classical [7,4,3] Hamming code, so one can correct at least up to 
2 erasures, which are the consequence of measurement failures. Thus, the encoded failure rate 
for measurements in the first level of encoding is 



assuming that classical processing needed for classical error correction is perfect. This yields 
a threshold of approximately 5 = 0.25, which guarantees the validity of both thresholds 
described above. 

These two cases, 5 = and 5 = ei,, are the extremal cases for the calculation of the threshold. 
In the limit where loss in the CSIGNs is much more likely than at the detectors, we have 6 = 
0, while in the limit where loss at the detectors is of the same order as loss in the CSIGNs, 
we have 6 = et. Because of the general construction of the CSIGN based on photon mode 
teleportation, there will be more than one detector involved, so the probability that there is a 
failure in the C S I GN teleportation will always be higher than the probability of a single detector 
failing, justifying this choice of upper and lower bounds. 




(5.26) 
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Lossy Error Model: Contour Plot of e/e/^^ as a function of e| and 5 




100 
50 
10 
5 
2 
1 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 



El 

Figure 5.5: Contour plot of -ffp as a function of and 6 under the lossy error model. The 
threshold is given by = j, and the greater this ratio is, the greater the gain for using erasure 
correction codes. 
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Figure 5.6: Plots of -jjj under the lossy error model, for extremal probabilities of detector failure, 
6 = and 6 = ej,. 



Chapter 6 
Simulation 



In order to verify the theoretical predictions about the recursion relation for the Steane code, a 
Monte Carlo simulation of the erasure correcting procedure was made following the principles 
previously used by Zalka OTl . which I review here for completeness. The results of the sim- 
ulations for the ideal and lossy error models are presented and contrasted with the theoretical 
predictions and calculations of Chapter|5l 

6.1 Background 

Fully simulating the states of the quantum computer, aside from having an exponentially high 
cost in resources, is superfluous for the determination of the performance of an error correcting 
code. Since we are only concerned with whether the error correction procedure is able to elim- 
inate all erasures, we can keep track of only that aspect of the computation. This is the same 
basic idea used by Zalka in [311, but applying various simplifications due to the simpler handling 
of erasures as well as new insights into fault-tolerant universal computation by various authors 

jHnmisi. 

While in the case of general errors one must keep track of how computation and error cor- 
rection propagates error operators that remain undetected, when dealing with erasures we simply 
need to keep track of where they appear during computation and erasure correction. The error 
operators that are kept track of are simply Pauli operators over the entire code block. 

It is crucial to note that this setup allows only for the simulation of gates that map the Pauli 
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group onto itself under conjugation, since, by assumption, all the erasure operators considered 
here are described by convex sums of Pauli operators - the Clifford group 62 described previ- 
ously. This is not a universal set for quantum computation, but it is enough for error correction. 
In order to make a universal set, one can include, for example, gates conditional on measurement 
outcomes. This was the method employed by Shor in the seminal paper on fault-tolerant compu- 
tation to demonstrate the construction of an encoded Toffoli for the Steane code ll26l . and it was 
further generalized to demonstrate the construction of other encoded gates such as the ^ gate and 
the controlled phase gate lfTni33l . 

Due to the fact that the encoded Toffoli can be built out of encoded Clifford gates, it is in- 
tuitive that a threshold for Clifford gates would also be a threshold for the Toffoli. This was 
analyzed and quantified carefully by Gottesman lIHl for the Steane code, where it was demon- 
strated that indeed the threshold for universal computation with Clifford gates plus the Toffoli is 
only ~ 5% lower than the threshold for Clifford gates. In the case of linear optics, the threshold 
would depend on the detector failure rate, since in order to construct these fault-tolerant gates 
necessary for universal computation measurements and operations conditioned on measurements 
are necessary. However, we have shown in Section l5.1.3l that the operation of the Clifford gates 
themselves, and the measurement of the syndromes already places enough constraint in the de- 
tector efficiency so that a result similar to the 5% approximation of the universal threshold given 
by Gottesman should still hold. 

6.2 Data Structures and Algorithms 

The simulation consists, in essence, of two arrays of bits of length 7, each position representing 
one of the 7 qubits that constitute a codeword in the Steane code. One of the arrays, called the 
X array, indicates the application of an X erasure on each qubit, while the other array, called the 
Z array, indicates the application of a Z erasure on each qubit - this is depicted in Figure 16.11 
and can be thought of as the side information mentioned in Section 12.11 One can represent full 
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Figure 6.1: Internal representation ofX®I®Z®E®I®I®Iin the Monte Carlo simulation. 

erasures as well, as X and Z commute, and 

Z(X(p)) = l(X(p)+ZX(p]Z) (6.1a) 

= ^ (p + XpX + ZpZ + YpY) (6.1b) 
= X(Z(p)) (6.1c) 
= E(p). (6.1d) 

In the simulation, a full erasure in qubit i is therefore represented by setting the ith bit in both 
the X and the Z arrays. 

As has been stated, the indication of any type of erasure over some qubit is by no means 
an indication that an error has indeed happened. The convex sum representation of the erasure 
superoperator can be interpreted, instead, as an indication that an error might have happened 
at that location. It is still necessary to measure the syndrome of the qubit erased to determine 
which of the possible Pauli operators has occurred. As described in Section l3^ the only sources 
of erasure that are considered here are the failures in teleportation of the CS IGN, and photon loss 
in the teleportation and detection of qubits. It has also been shown that these are all erasures, so it 
is always known when the data is corrupted by these failures, and one may simply try to generate 
a large number of ancillae in the desired state and discard all the ones that contain erasures. 

In the simulation, we assume the erasure correction procedure is applied after a single appli- 
cation of an encoded CSIGN gate resulting in erasures. Because the CNOT and the Hadamard 
are transverse on the Steane code (see Section 14.1.1 1) , so is the CSIGN, and it realizes the en- 
coded CSIGN itself. In that case, we assume each of the gates fails independently at the target 
and at the control. Under the ideal error model, the only type of failure that may occur is a 
Z measurement, and it is localized to either the control or the target independently. Under the 
lossy model, however, a teleportation failure at the CSIGN induces a full erasure locally and a 
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Z erasure at the gate counterpart - for example, a failure at the target causes a full erasure at the 
target and a Z erasure at the control. Thus, the error pattern on a code block depends not only 
on the failures that occurred locally, but also on the failures that occurred on a block with which 
it was interacting, but again, because of the rule given by (16.11) . this is easy to compute. Under 
both erasure models one simulates the initial state of the block by flipping a biased coin for each 
qubit to determine whether it should have an erasure or not, and in the case of the lossy error 
model, there are two passes: one determines whether there should be a full erasure, and the other 
determines whether it should have a Z erasure if it did not have a full erasure. 

Once the erasure pattern has been generated, multiple rounds of single erasure correction are 
attempted until either the code block is erasure free, or the erasure pattern is uncorrectable. Since 
each single erasure correction step requires the measurements of a weight four stabilizer operator 
(in the case of the Steane code), any of these gates may fail at either the control or the target side, 
independently, so once again biased coins are flipped. In this case, the difference between the 
two error models is a little more salient. For the ideal model, we can simply count how many 
failures there were on the ancilla side of the stabilizer measurement in order to determine if the 
error correction will succeed or fail, and then determine which gates failed on the data side so 
that the erasure pattern may be updated for the next correction attempt. Under the lossy model, 
however, one must keep track of which CSIGNs with the ancilla failed, because those failures 
will induce erasures on the data side, which must be remembered for the next correction attempt, 
as well as keeping track of which failures occurred on the data side. Which qubits will need to 
interact with the ancilla is determined by the erasure pattern, as described in Sections 14.1.21 and 

Km 

The lossy model has another distinction in that we allow for erasures to occur at the measure- 
ment of the ancilla, although at a rate 5 different from the gate failure rate. In reality ei and 6 
are not independent, since constructions of the C S I GN gate itself depend on the measurement of 
ancillary photonic modes. The exact relationship is highly dependent on which CS IGN construc- 
tion is being used, and since each has its own advantages, it was decided to make the simulation 
take ei and 5 as independent parameters that can be adjusted as desired. 

For each gate failure rate the simulation is run until at least one thousand blocks are found to 
be uncorrectable (for the low probabilities of failure), or until at least a large number of blocks, 
proportional to -, have been simulated (for high probabilities of failure). Statistics, such as 
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Ideal Error Model: log-log scale e/e;'^' vs. e. 
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Figure 6.2: Probability of failure for teleportation failure rate e using 6 erasure correction at- 
tempts. 

how many blocks were simulated, how many were uncorrectable, etc. are output into a file, and 
the process is repeated for a different random number generator seed. After 20 such iterations, 
the individual statistics are compiled to obtain a mean and a standard deviation for the encoded 
failure rate at every gate failure rate (and detector failure rate). 



6.3 Results 

6.3.1 Ideal error model 

As can be observed in Figure the simulation matches the predicted values very closely. Be- 
cause of repeated trials with different seeds for the random number generators, the error bars for 
the simulation results are very small, and we can take the variance in the results to be insignifi- 
cant. 
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Lossy Error Model: log-log scale E|/E|'^' vs. E| 
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Figure 6.3: Probability of failure given 5 = with 20 erasure correction attempts. 

6.3.2 Lossy error model 

Again, as can be seen in Figures lOl and the error bars in the simulation are so small as to 
be insignificant, but in this case there is a discrepancy between the simulation results and the 
calculated values for probabilities below the threshold. 

There are a number of possible reasons for this behavior, but the main source of error appears 
to be the recursion relation itself and using computers to numerically approximate its value at 
different gate failure probabilities. While it is necessary to perform the 20 correction attempts in 
order to have the same precision in the leading terms of the recursion relation as in the ideal case, 
that makes the recursion relation extremely long, with very large integer coefficients, both posi- 
tive and negative. By construction, this recursion relation must stay bounded to values between 
and 1 . However, while the mathematical packages have infinite precision for integer handling. 
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Figure 6.4: Probability of failure given 6 = ei, with 20 erasure correction attempts. 



the same is not true for floating point calculations, necessary in the evaluation of these functions. 
The finite precision of these calculations used to generate these plots should potentially be the 
main source of error. The recursion relations themselves are correct. 

A different approach to verify the correctness of the recursion relations would be to fit poly- 
nomial curves to the simulation results. However, in order to do so reliably, one needs to run 
simulations at extremely low probabilities, increasing the error on the curves significantly sim- 
ply because not enough statistical sampling is done. This is the same problem observed by Steane 
and it is a common problem of Monte Carlo simulations at low event probabilities. 



Chapter 7 
Conclusion 



Through theoretical calculations and verification by simulation, we have found that in an error 
model consisting only of ideal teleportation failures, there is a threshold of ^ 0.1 15 = 11 .5%, 
and that for an ideal model consisting of photon loss and replacement during teleportation, the 
threshold is at least ei, ^ 0.01 78 = 1 .78% as long as the detector failure rate at most equal to the 
gate failure rate. A mixed error model consisting of both types of failure would have a threshold 
somewhere in between these two as long as the total probability of failure e is bounded by these 
two values, that is 

eiKcKCi. (7.1) 

In the limit where the photon loss dominates, one expects the threshold to be closer to ei^, while 
in the case where ideal failures dominate, the threshold would be higher and closer to e|. 

The results presented here need to be considered carefully. The correction procedure used 
for the threshold calculation can correct either a single general error at an unknown location, or 
two (and at most three) erasure patterns. It cannot handle a mixed error model involving both 
errors and erasures. For realistic quantum computation we need to consider all of these sources 
of error, so one would need to consider a modified correction procedure. 

However, the main result to be taken from this work is that the recursion relation for the 
erasure correction procedure can be derived exactly, and this may lead to new techniques of 
analyzing the performance of erasure correcting codes, or even general error correcting codes. 
The erasure correcting procedure used here is not fully optimized, since it was tailored to simplify 
the exact analysis, but there are many avenues that may be taken to improve these bounds on the 
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erasure threshold, including different syndrome measurement techniques lESlI^ . It is important 
to note that the thresholds presented here are well above the general error thresholds presented in 
ll29ll3n . and although the threshold for the ideal case if substantially worse than the eum = 0-5 
threshold given in ifTTll . my result considers a more general and more realistic error model. 

7.1 Future Work 

In order to make the results presented here of direct applicability to linear optics quantum com- 
puting, one will have to take the different construction of quantum gates into account, and de- 
termine how ei depends on the detector efficiency 5. There has been some previous work in 
that area ll20l . but none that I am aware relating detector efficiency directly to a threshold for 
quantum computing. This relationship would be of extreme importance in determining what the 
requirements are for the performance of single photon detectors. The results presented here are 
the first step in that direction. 

As mentioned before, one would need to consider larger codes or different error correction 
procedures in order to be able to correct both erasures and general errors. Recently a large survey 
of the performance of different CSS codes has been published ll29ll . and it has even been pointed 
out how the [[23, 1 , 7]] Golay code provides a better threshold for general errors than the Steane 
code. A modification of the simulation used to obtain the results presented here should be straight 
forward, especially with the insights into how to perform erasure correction by using information 
about the automorphism of the code. There is no comprehensive study of how different codes 
compare for an erasure error model, although there have been many papers dealing with erasure 
correction codes ifT^ [1311 . so a survey of erasure thresholds along the lines of could be 
elucidating. 

Although for the error model presented here it would be unreasonable to choose anything but 
a CSS code, for other error models that may have erasures at any gate the same may very well 
not be true. A modification of the simulation could be done based on ideas on how to encode 
stabilizer codes in binary string (presented in ll25ll ). but how to modify the code to make it general 
and flexible enough for different error models is more involved. A natural first step would be to 
simply use general error models for gate and memory error so comparisons can be made between 
known stabilizer codes and the CSS codes presented in lE^ . 
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Finally, the most interesting aspect of this research, in my opinion, has been the calculation 
of the recursion relations by describing the erasure correction procedure as a Markov chain. 
An investigation into how to extend this work into the general error model could yield precise 
insights into what governs the relationship between different code parameters, such as the weight 
distributions, the automorphisms of the code, etc. and the thresholds obtained with that code. 



Appendix A 

Calculation for Z measurement Threshold 



Following the procedure described in Section |231 we take ei to be the probability of a Z mea- 
surement occurring due to teleportation failure. 

The error free pattern, that is, a measurement pattern of weight 0, is an absorbing state, so the 
only non-zero probability of transition is into itself - nothing needs to be corrected. The same is 
true for the uncorrectable measurement patterns, represented by the fail state. 

Pr(i|0) = 5,,o (A.l) 
Pr(iM) = (A.2) 
Fr [faillfail] = 1. (A.3) 

Regardless of what the current state, the probability of correcting a single Z measurement is 
given by the probability that there are no failures at all is the correction procedure, that is, all 
three teleportations succeed, or 

Pr[success) = (1 - ei)^. (A.4) 

The probability of staying in an measurement pattern of the same weight is since any tele- 
portation failure in the correction procedure prevents the restoration of the previous measured 
qubit. 

Many of the transition probabilities in the measurement correction procedure used here are 
zero. Exactly which transition are possible given the procedure follows from some simple rules. 
Given some Z measurement pattern of weight i 
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• only one measurement is corrected at a time, resulting in an measurement pattern of weight 
i — 1 in the best case. 

• otherwise, the weight of the measurement pattern can only increase or stay the same (Mur- 
phy's Law for Z measurement correction). 

Using general cases described here, we calculate all the non-zero probabilities of transition, 
omitting detailed description of how they were obtained for brevity. All probabilities of transition 
equaling zero are omitted. 



A.l Transitions from weight 1 patterns 



Pr(0|l) 



Pr[success] 



(A.5a) 




(A.5b) 



(A.5d) 



(A.5c) 



A.2 Transitions from weight 2 patterns 



Pr(1|2) 



Pr[success) 



(A.6a) 




(A.6b) 



(A.6c) 



A.3 Transitions from correctable weight 3 patterns 



Pr(2|3) 
Vr{fail\3) 



= Pr[success) 



(A.7a) 
(A.7b) 
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A.4 Initial distribution 

Since we are only considering a single type of error - namely, unintentional Z measurements due 
to teleportation failure - this calculation is quite simple, yielding 

Pr(j) = ^.]4{^-e,V-\]<3 (A.8a) 



Pr(3] = ^(sViH-et)" (A.8b) 



Fr {fail) - 1(^^40 + (A.8c) 



Appendix B 

Calculation for Full Erasure Threshold 



Following the procedure described in Section l231 we take ei^ to be the probability of a full erasure 
occurring, and 6 to be the probability that a single detector will fail. The general idea is to adapt 
the procedure used for correcting Z measurements and use it to correct Z erasures, and use the 
standard fault-tolerant stabilizer measurement of Section l2.6.1l to correct X erasures. 

The probability that there will be a failure in any of the detectors in the measurement of the 
ancilla (in the fault- tolerant stabilizer measurement) is 

Vxiancilla det.] = ^ -[^ - 5)^, (B.l) 

since all Z stabilizer operators of the Steane code have weight 4. The error free state [0, 0] is an 
absorbing state, so the only non-zero probability of transition is into itself. The same is true for 
the uncorrectable errors, represented by the fail state. 

Pr([iJ] I [0,0]) = 5t,o5j,o (B.2) 
Pr([iJ]M) = (B.3) 
Fr {fail\faU) = 1. (B.4) 

Regardless of what the current state is, the probability of correcting the single erasure we 
were targeting is given by the probability that there are no failures at all. In the case of an X 
erasure, that means 

VrisuccesslK) = (1 — e;)^(1 — Fr[ancilla det.]] . (B.5) 
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In the case of a Z erasure, we need to consider the qubit teleportation described in Figure 14.41 
There we have a single entangling gate, the 90° rotation about Z ® Y, that like the CSIGN 
is implemented through mode teleportation, so a photon is lost with the same probability £]_. 
Then each of the two qubits is measured with an independent probability of failure 5. Thus, the 
probability that a photon is lost in this teleportation is given by 



We also need to measure the Z erasure in order to apply the same correction procedure as in the 
Z measurement, yielding 



When attempting to correct a Z, the probability of staying in the same state is 0, since any 
failure increases the weight of the erasure pattern and re-randomizes the phase of the qubit with 
a Z erasure, while failure during the initial measurement transforms a Z into a E. On the other 
hand, when there are only full erasures and an X must be corrected (recall X are always corrected 
first, by design), in order to stay in the same erasure pattern, either only the measurement must 
fail, the CS IGN on the target failed on the data side, or the CS IGN failed only on the ancilla side 



Many of the transition probabilities in the erasure correcting procedure used here are zero. 
Exactly which transitions occur with non-zero probability, given the procedure, follows from 
some simple rules. Given some erasure pattern described by [m, n] where m is the number of 
full erasures and n is the number of Z erasures, 

• the type of erasure (X or Z) that is most abundant is corrected first. 

• only one erasure is corrected at a time, but when a full erasure is corrected, it becomes a Z 
erasure, resulting in [m — 1 , n + 1] in the best case. If a Z erasure is being corrected, the 
result is [m, n — 1] in the best case. 




(B.6) 



Fr{success\Z] = (1 -5) (1 -Vr{loss]f . 



(B.7) 



Pr(stay E) = (1 - eO^[(1 - eif Pr [ancilla det.) + (1 - et]ej 



(B.8) 



• a Z erasure may be transformed into a full erasure if, when the qubits is measured in the Z 
basis, the detector fails. 
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• failure during a Z erasure correction results in the Z erasure staying as a Z erasure, but 
introducing full erasures in the other three qubits in the subsystem. Thus, only the number 
of full erasures can increase due to failure during a Z erasure correction - by a minimum 
of 1 and a maximum of 3. 

• Failures during an X erasure correction results in the number of either types of erasures 
increasing or stay the same. 

Using the general cases described here, we calculate all the non-zero probabilities of tran- 
sition, omitting detailed description of how they were obtained for brevity. All probabilities of 
transition equaling zero are omitted, and the probabilities of transition to the fail state are chosen 
such that the probabilities of transition from any given state sum up to 1 . 



B.l Transitions from [0,1] 



Pr([0,0]|[0,1]) 
Pr([1,0]|[0,l]) 



Pr(success|Z) 



5 



(B.9a) 
(B.9b) 




(B.9d) 



(B.9c) 
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Pr([0,l] 
Pr([1,0] 

Pr([1,l] 



Pr([2,0] 
Pr([2 
Pr([1,2] 
Pr([3,0] 



1]|[1 



[1,0]) 
[1,0]) 

[1,0]) 



[1,0]) 
,0]) 
[1,0]) 
[1,0]) 



Pr( success |E) 
Pr(stay E) 

(l-eO'(^]e,(1-eO' 



^0e.(1-eO^0e.(l-eO 

^(l-eO^Q)^f(l-eO 
2/3 



3 V2 



ef(l-ei) 



(B.lOa) 
(B.lOb) 

(B.lOc) 

(B.lOd) 

(B.lOe) 

(B.lOf) 

(B.lOg) 
(B.lOh) 



B.3 Transitions from [1,1] 



Pr([1,0]|[1,1]) 
Pr([2,0]|[l,1]) 

Pr([2,1]|[l,l]) 



Pr(success|Z) 
5 

(1 -S)^K j Pr(loss)[l -Pr(loss)]^ 



B.4 Transitions from [2,0] 



(B.lla) 
(B.Ub) 

(B.llc) 



Pr([l,1]|[2,0]) 
Pr([2,0]|[2,0]) 

Pr([2,l]|[2,0]) 
Pr([3,0]|[2,0]) 



Pr(success|E) 
Pr(stay E) 



2/3 



3 VI 



ei(l -eO"(l -eO^ 



(B.12a) 
(B.12b) 

(B.12c) 
(B.12d) 
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B.5 Transitions from [0,2] 

Pr([0,1]|[0,2]) = Pr(successlZ) (B.13a) 
Pr([l,1]|[0,2]) = 5 (B.13b) 

Pr([l,2]|[0,2]) = (1 -6)^(^^^ Pr(loss)[l -Pr(loss)]2 (B.13c) 

B.6 Transitions from [3,0] 

Pr([2,1]|[3,0]) = Pr(successlE) (B.14a) 

Pr([3,0]|[3,0]) = Pr(stayE) (B.14b) 

B.7 Transitions from [0,3] 

Pr([0,2]|[0,3]) = Pr(successlZ) (B.15a) 

Pr([l,2]|[0,3]) = 6 (B.15b) 

B.8 Transitions from [2,1] 

Pr([2,0]|[2,1]) = Pr(successlZ) (B.16a) 

Pr([3,0]|[2,1]) = 5 (B.16b) 

B.9 Transitions from [1,2] 

Pr([l,1]|[l,2]) = Pr(successlZ) (B.lVa) 

Pr([2,1]|[1,2]) = 5 (B.lVb) 
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B.IO Initial distribution 



The initial distribution can be computed easily by noting the simple facts that 



Z(E(p)) 
E(Z(p)) 
Z(Z(p)) 
E(E(p)) 



E(p) 
E(p) 
Z(p) 
E(p). 



(B.lSa) 
(B.lSb) 
(B.lSc) 
(B.lSd) 



Assuming that each C S I GN has two different failure modes - control and target failure - we 
first calculate the probability of getting full erasures for the desired pattern. Once that has been 
calculated, we calculate the probability of getting Z erasures on the remaining unaffected qubits. 
Following this procedure, one finds 



Pr( 
Pr( 

Pr( 

Pr( 

Pr( 

Pr( 

Pr( 
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Pr( 

Pr( 



0,0] 
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(1-eO 



14 



e^(1 -et)^(1 -et)^ 



4 
5 
4 
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4 
5 
4 
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{^-e,v(l]e![^-e,r 



(B.19a) 
(B.19b) 

(B.19c) 

(B.19d) 

(B.19e) 

(B.19f) 

(B.19g) 

(B.19h) 

(B.19i) 

(B.19j) 
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and the probability of seeing an uncorrectable erasure is simply one minus the sum of these 
probabilities. 
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